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ABSTRACT 

> 

Laminar  fluid  flow  in  a simplified  fluidic  amplifier  is  analyzed  using 

i 

finite  difference  numerical  techniques.  Numerical  solutions  to  the  two  dimen- 
sional incompressible  Navier-Stokes  equations  are  obtained  in  terms  of  stream 
function  and  vorticity.  These  solutions  are  presented  as  contour  plots  of 
lines  of  constant  stream  function.  The  effect  of  various  flow  parameters,  such 
as  Reynolds  number,  outlet  port  loading,  and  deflection  of  the  power  jet  are 

* 

examined.  Finite  difference  solutions  are  also  obtained  for  the  pressure  dis- 
tribution. The  problems  associated  with  obtaining  accurate  numerical  pressure 
solutions  are  discussed  in  detail.  In  particular,  the  problems  associated  with 
an  inconsistent  formulation  of  this  Neumann  problem  are  described.  A large- 
scale  laminar  flow  amplifier  was  designed  and  constructed  so  that  the  geometry 
1 ‘ matched  that  of  the  numerical  simulation.  Dye  was  injected  into  water  flowing 

in  the  amplifier  and  streamline  photographs  were  obtained.  Good  agreement 
existed  between  the  photographs  and  the  numerically  predicted  streamlines. 

Extensions  of  the  model  aimed  at  producing  a working  tool  for  use  in  design 
of  fluidic  amplifiers  are  suggested. 
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EXECUTIVE  SUMMARY 

The  basic  components  of  a beam  deflection  fluid  jet  amplifier  are  a 

power  jet  which  issues  from  a nozzle  and  impinges  on  symmetrically  placed 

recovery  passages,  and  control  passages  which  act  on  the  side  of  the  power 

jet  so  as  to  give  the  resulting  stream  a lateral  deflection  relative  to 

★ 

the  line  of  symmetry.  The  recent  Fluidic  State  of  the  Art  Symposium  pro- 
vides excellent  background  for  the  reader  who  is  unfamiliar  with  the  field. 

Static  Output  Characteristics 

Resulting  from  early  work  under  this  contract  it  was  found  that  when 
a power  jet  impinges  on  a single  receiver  passage,  the  recovery  pressure 
and  flow  at  the  output  port  of  the  receiver  could  be  calculated  only  when 
receiver  loading  allowed  significant  through  flow  of  fluid.  A method  in- 
volving the  use  of  moving  "cowl -streaml ines"  with  a lumped  control  volume 
model  assuming  inviscid  flow  in  the  control  volume  region  was  proposed  as 
an  improved  and  viable  way  to  predict  output  performance. 

This  work  was  later  extended  to  the  case  of  a jet  interacting  with 
multiple  receivers  and  by  careful  selection  of  the  location  of  the  cowl- 
streamlines,  again  it  was  possible  to  predict  the  output  recovery  in  the 
receiver  system.  Although  a wide  range  of  jet  flow  Reynolds  Numbers  and 
receiver  geometries  were  not  investigated,  the  predictive  power  of  the 
model  was  adequate  (less  than  10%  error)  for  the  cases  considered. 

Amplifier  Jet-Vent  Dynamics 

Moving  to  the  dynamics  of  a complete  amplifier,  early  studies  had 
shown  that  the  bandwidth  of  a beam-deflection  amplifier  was  limited  by  the 

* Fluidic  State  of  the  Art  Symposium--F!arry  Diamond  Laboratories,  September 
1974. 
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inertia  of  the  side  vent  fluid.  One  of  the  unknown  quantities  in  the 
development  of  theoretical  models  for  predicting  bandwidth  was  the  coefficient 
relating  to  the  division  of  the  power  jet  in  the  side  vent  areas  to  main  jet 
deflection. 

This  coefficient  depends  on  output  loading  and  receiver  configuration. 

Thus  utilizing  the  work  of  the  early  studies  on  recovery  and  spillover 
characteristics,  a lumped  system  dynamic  model  including  jet  transport  delay 
effects  was  formulated.  The  results  of  the  work  showed  an  almost  quadratic 
reduction  in  bandwidth  as  the  side  vent  lengths  were  increased  but  only 
compared  accurately  with  experimentally  measured  bandwidths  for  short  side 
vents.  Higher  frequency  resonances  were  seen  and  were  attributed  to  jet 
lateral  dynamic  modes  in  which  the  apparent  signal  transmission  speed  in 
a turbulent  jet  was  found  to  travel  at  nearly  one-fourth  of  the  power  jet 
nozzle  flow  velocity. 

Numerical  Models  for  Laminar  Interacti on  Regions 

Finally,  to  investigate  the  feasibility  of  using  new  efficient  numerical 
codes  for  the  solution  of  the  laminar  flow  field  equations,  a study  was  under- 
taken of  a non-vented  amplifier  without  control  jet  action.  The  equations  of 
laminar,  two  dimensional,  incompressible  flow  were  solved  using  finite  difference 
methods.  The  solutions  for  distribution  of  stream  function  values  showed  good 
agreement  with  experimentally  observed  flow  patterns. 

Numerical  methods  were  compared  for  calculating  internal  pressure  dis- 
tributions with  asymmetric  loading  at  the  receiver  ports  so  that  output  port 
pressure-flow  characteristics  could  be  predicted.  The  accuracy  of  the  pressure 
distribution  simulation  was  not  satisfactory  and  indicated  that  problems  can 


arise  because  of  the  need  to  solve  a Poisson  type  differential  equation  with 
Neumann  boundary  conditions.  An  inconsistency  in  the  equations  can  be  found 
under  some  conditions  and  numerical  errors  can  reduce  the  final  accuracy  in 
the  predicted  pressure-flow  relationships  at  the  amplifier  ports. 

The  details  of  the  numerical  modeling  research  are  presented  in  this 
final  technical  report. 
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CHAPTER  I 

INTRODUCTION 

The  proportional  amplifier  is  one  of  the  basic  components  of 
many  fluidic  systems.  Since  the  early  1960's  a great  deal  of  work 
has  been  done  modifying  and  optimizing  the  design  of  such  amplifiers, 
understanding  the  fluid  mechanics  involved,  and  developing  techniques 
for  the  use  of  amplifiers  in  control  circuits. 

As  with  most  engineering  components  fluidic  devices  have  both 
advantages  and  disadvantages  [10,16].  One  important  advantage  is  the 
ability  of  properly  designed  fluidic  sys tents  to  function  in  harsh 
environments,  especially  in  comparison  to  electronic  systems  designed 
to  perform  similar  functions.  Often  having  no  moving  mechanical 
parts,  fluidic  devices  are  resistant  to  failure  due  to  mechanical 
shock  or  vibration.  They  can  be  designed  with  low  cost  to  withstand 
high  temperatures  or  high  radiation  fields. 

Unfortunately,  there  are  attendant  disadvantages.  Relatively 
low  speed  of  operation  compared  to  electronics  results  from  both  the 
relatively  slow  acoustic  propagation  speed  in  the  fluid  and  the  time 
required  for  fluid  movement.  The  tendency  of  fluid  flows  to  become 
unstable  generates  turbulent  fields.  This  turbulence  gives  rise  to 
high  noise  content  in  the  sensing  and  control  signals.  Also,  ampli- 
fier gains  and  input  impedances  are  low  compared  to  corresponding 
electronic  devices. 
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Perhaps  a more  subtle  problem  facing  the  development  of 
fluidics  is  the  lack  of  a completely  adequate  model  for  predicting 
the  performance  of  components.  The  characteristics  and  usefulness 
of  all  fluidic  devices  depend  on  the  interaction  of  one  or  more 
fluid  jets.  Most  modeling  techniques  to  date  have  relied  on  assump- 
tions of  flow  conditions  for  these  confined  jets  and  have  met  with 
success  only  in  specific  cases.  The  search  for  a more  general  way 
of  modelinq  the  internal  fluid  flow  interactions  in  fluid  amplifiers 
has  been  the  motivation  for  tiris  work. 

A.  Beam  Deflection  Proportional  Amplifiers  and  Characteristic  Curves 

Consider  in  particular  a beam  deflection  proportional  ampli- 
fier. Although  the  fluid  mechanics  involved  in  such  a device  are 
complex,  the  basic  idea  is  very  simple.  A typical  amplifier  is 
shown  in  Fig.  1.1.  The  differential  signal  which  is  to  be  amplified, 
either  a pressure  or  a flow,  is  established  as  an  input  between  the 
control  ports. 

A high  energy  supply  jet  (typically  having  ten  times  the 
stagnation  pressure  level  of  the  control  flows)  enters  the  amplifier 
through  the  power  nozzle.  It  is  deflected  by  the  differential  pres- 
sure and  momentum  of  the  control  flows.  This  main  flow  continues 
across  the  amplifier  cavity  at  some  angle  to  the  cavity  centerline 
and  impinges  on  the  receivers.  The  vents,  receivers,  and  center 
dump  collect  varying  amounts  of  flow  and  energy  depending  on  the 
degree  of  deflection  of  the  power  jet.  Thus,  a differential 
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pressure  is  produced  at  the  receiver  outlets.  This  differential  is 
directly  related  to  the  control  pressure  differential  but  of  larger 
magnitude;  that  is,  the  original  differential  signal  has  been  ampli- 
fied. 

A large  variety  of  fluid  flow  phenomena  are  involved  in  this 
process.  Consider  the  main  jet  as  it  flows  across  the  amplifier.  To 
begin  with,  the  jet  flow  can  be  either  laminar  or  turbulent.  As  the 
jet  leaves  the  power  nozzle  a velocity  profile  begins  to  develop  due 
to  the  viscous  or  turbulent  shear  on  the  flow.  Quickly  this  profile 
is  modified  by  interaction  with  the  control  flows.  Further,  the 
geometry  in  the  control  region  will  effect  the  developing  profile. 

Of  particular  interest  are  tne  positioning  of  the  control  nozzles 
and  the  amplifier  aspect  ratio,  which  is  the  ratio  of  the  power 
nozzle  width  to  the  amplifier  depth.  Once  in  the  cavity  region  the 
jet  is  acted  on  by  shear  forces  and  by  the  recirculation  flow  from 
the  amplifier  vents.  As  it  impinges  on  the  receivers  the  jet  profile 
is  destroyed  by  stagnation  against  the  wall  arid  by  division  of  the 
flow  between  the  receivers,  vents,  and  center  dump.  This  flow  divi- 
sion is  extremely  important  but  little  understood. 

The  staging  and  design  of  beam  deflection  amplifier  circuits 
are  described  in  detail  by  Belsterling  [3j.  This  procedure  will  not 
be  described  here  except  to  point  out  that  it  is  necessary  to  know 
the  "characteristic"  curves  of  a given  amplifier.  These  curves 
graphically  describe  how  the  amplifier  will  respond  to  various 
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> control  pressure  differentials  and  what  the  relationship  is  between 

both  the  control  flows  and  pressures  and  the  output  flows  and  pres- 
sures. 

For  example,  a typical  static  input  characteristic  curve  for 
one  control  port  of  an  amplifier  is  shown  in  Fig.  1.2.  A similar 
curve  would  exist  for  the  other  port.  The  curve  relates  the  flow 
< through  the  control  port  to  the  pressure  at  the  port. 

The  influence  of  the  control  pressure  differential  on  the 
receiver  output  pressures  is  commonly  shown  in  one  of  two  ways.  The 
output  pressure  for  one  of  the  receivers  can  be  plotted  as  a function 
of  the  control  pressure  differential,  as  in  Fig.  1.3.  Alternatively, 
the  pressure  differential  at  the  receivers  can  be  shown  as  a func- 
tion of  the  control  pressure  differential,  as  in  Fig.  1.4. 

The  output  characteristics  for  the  amplifier  are  typically 
plotted  as  shown  in  Fig.  1.5.  Here  the  flow  in  one  receiver  is 
plotted  as  a function  of  the  output  pressure  of  that  receiver.  Sep- 
arate curves  are  given  for  various  values  of  the  control  pressure 
differential.  All  of  the  above  characteristic  curves  are  obtained 
experimentally  as  explained  by  Belsterling  [3]. 

While  the  particular  amplifier  shown  in  Fig.  1.1  is  vented, 
other  amplifiers  are  also  used  in  which  the  side  vents  are  closed. 

This  eliminates  the  need  for  a separate  drain  to  handle  the  vent 
flows.  Naturally,  non-vented  amplifiers  are  only  used  when  the 
loading  at  the  receivers  allows  continuous  through  flow.  Otherwise 
the  flow  in  the  power  jet  would  be  stopped. 
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Figure  1.2.  Typical  Input  Characteristic  Curve 
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Figure  1.3.  Typical  Single-sided  Pressure 
Gain  Characteristic  Curve 
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Optimization  of  design  and  performance  characteristics  has 
been  accomplished  in  the  past  by  trial  and  error  experimental  proc- 
esses coupled  with  some  assistance  from  limited  theoretical  models. 
These  models  have  employed  control  volume  analyses  [11].  The  major 
drawback  with  these  models  is  the  assumptions  which  must  be  made  as 
to  velocity  and  pressure  distributions  within  the  internal  field  of 
flow. 

The  aim  of  this  study  was  to  examine  the  problems  and  inves- 
tigate the  feasibility  of  using  finite  difference  techniques  for 
predicting  the  internal  flow  and  performance  characteristics  of  a 
simplified  beam  deflection  amplifier.  With  these  techniques  many  of 
the  assumptions  of  the  control  volume  analyses  can  be  avoided. 

The  choice  of  the  particular  amplifier  to  be  modeled  was 
governed  by  several  considerations.  Because  of  the  improved  signal 
to  noise  ratios  in  laminar  amplifiers  [13],  a laminar  flow  device 
was  chosen  for  the  study.  Also,  in  order  to  keep  the  scope  of  the 
investigation  within  bounds,  a simplified  device  was  chosen  which 
resembles  a non-vented  amplifier  having  no  control  passages.  The 
reason  for  eliminating  control  passage  effects  at  this  stage  was  to 
concentrate  on  the  interactions  between  the  power-jet,  the  receivers, 
and  the  side  circulation  regions  of  the  device.  The  effects  of  jet 
deflection  were  modeled  here  by  a simple  angular  deflection  of  the 
power  jet  relative  to  the  device  centerline.  It  was  felt  that 
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control  flow  effects  could  be  included  later.  Finally,  a high  as- 
pect ratio  has  been  chosen  for  the  amplifier  in  this  investigation. 
While  this  may  not  be  the  best  situation  practically  for  fluid  ampli- 
fiers, it  does  allow  the  assumption  for  now  of  a two  dimensional  flow 
field.  Extensions  of  the  work  to  three  dimensional  flow  could  be 
made  in  principle. 

C.  Method  of  Approach 

The  study  divides  into  three  major  parts.  First,  the  stream 
function  and  vorticity  transport  equations  are  solved  using  finite 
difference  techniques.  Basically  the  model  assumes  that  the  ampli- 
fier flow  is  a laminar,  two  dimensional,  incompressible,  viscous 
flow.  Second,  based  on  the  calculated  flows  finite  difference  solu- 
tions are  obtained  for  the  pressure  distribution  in  the  amplifier. 
Finally,  experimental  flow  visualizations  are  obtained  through  dye 
injection  into  water  flowing  in  a large  scale  laminar  amplifier 
model.  These  visualizations  provide  qualitative  experimental  com- 
parisons for  evaluation  of  the  numerically  predicted  flow  patterns. 


CHAPTER  II 


DEVELOPMENT  OF  THE  NUMERICAL  MODEL 

This  chapter  begins  with  presentation  of  the  governing  set  of 
partial  differential  equations  for  the  fluid  flow  being  considered. 
These  continuum  equations  divide  naturally  into  two  separate,  but  re- 
lated, systems  of  finite  difference  equations;  i.e.  the  flow  field 
equations  and  the  pressure  field  equations.  Since  analytical  solu- 
tions for  these  problems  are  not  feasible,  numerical  solutions  using 
finite  difference  techniques  were  sought.  Appropriate  differencing, 
boundary  conditions,  method  of  solution,  and  convergence  criterion 
will  be  presented  for  both  systems.  In  addition,  for  the  pressure 
equation  the  problem  of  inconsistency  will  be  examined. 

A.  Continuum  Equations 

The  physical  process  being  modeled  is  that  of  a viscous  fluid 
moving  through  a confined  flow  region.  Since  Mach  numbers  are  small, 
the  flow  can  be  assumed  incompressible  and  with  the  restriction  to 
suitably  low  flow  velocities  the  assumption  of  laminar  flow  in  all 
regions  may  be  made.  Further,  it  is  assumed  that  the  fluid  proper- 
ties are  invariant  and  that  body  forces  are  negligible.  Finally,  it 
is  assumed  that  a two  dimensional  model  will  give  a reasonable  repre- 
sentation of  the  physical  process. 
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Under  these  assumptions  the  motion  of  the  fluid  is  described 
by  the  following  form  of  the  Navier-Stokes  and  continuity  equations 
[25].  In  terms  of  the  Cartesian  coordinates  x arid  y and  time  t, 
these  equations  are 
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where  u and  v are  the  x and  y velocities,  respectively,  P is  the 
pressure,  and  p and  v are  the  fluid  properties  density  and  kinematic 
viscosity,  respectively. 

In  order  to  simplify  Eqs.(2.1),  (2.2),  and  (2.3)  for  our  pur- 
poses here,  it  is  convenient  to  make  a change  of  variables  and  to 
eliminate  the  pressure  from  the  resulting  set  of  equations.  In 
addition,  throughout  this  work  the  variables  will  be  nondimension- 

alized  using  the  average  inflow  velocity  U,  the  width  b of  the  inlet, 

2 

and/or  the  quantity  pU  (see  Table  2.1).  The  variable  transforma- 
tion introduces  the  variables  stream  function  iJj  and  vorticity  r, 
which  are  defined  by  the  relations 
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Table  2.1.  Non-Dimensional  Variables 

Dimensional  Dimensionless 

Variable  Ratio 


A 

U 

u 

= u/U 

V 

V 

= v/U 

X 

X 

= x/b 

A 

y 

y 

JD 

< >> 

II 

p 

p 

= P/pU2 

? 

r 

S 

= S/(U/b) 

= i/(ub) 

Finally,  the  pressure  is  eliminated  by  differentiating  Eq.(2.1)  with 
respect  to  y and  Eq.(2.2)  with  respect  to  x.  Upon  subtraction  of  the 
two  resulting  equations  the  derivatives  of  pressure  cancel  and  the 
vorticity  transport  equation  is  obtained.  One  can  rearrange  this 
equation  into  conservative  form  [20]  and  obtain 


K - 3(uc)  a(vd  + J_  ^ 

« - ' 3x  ay  Re  [^7  ^7 

(advection  terms)  (diffusion  terms) 
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(2.6) 


(2.7) 


where  the  Reynolds  number.  Re,  is  Ub/v. 


With  application  of  suitable  boundary  conditions  this  system 
of  coupled  partial  differential  equations  defines  the  flow  field. 
Solution  for  the  flow  field  using  these  equations  is  simplified  com- 
pared to  use  of  the  original  set  because  the  pressure  terms  have 
been  eliminated. 

When  a pressure  field  solution  is  desired,  't  must  be  recon- 
structed from  the  flow  field  solution.  Pressure  distributions  can  be 
obtained  by  solution  of  a Poisson  type  partial  differential  equation 
which  is  also  derived  from  Eq . (2.1)  and  Eq . (2.2)  [20].  Differenti- 
ating Eq.(2.1)  with  respect  to  x,  differentiating  Eq . (2.2)  with 
respect  to  y,  adding,  eliminating  terms  through  use  of  continuity, 
and  transforming  variables  results  in  the  elliptic  equation 
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It  should  be  emphasized  that,  although  the  equations  are  separate, 
the  pressure  field  equation  is  dependent  on  having  a solution  for 
the  flow  field.  This  dependency  is  reflected  both  in  the  derivatives 
on  the  right-hand  side  of  Eq . (2.8)  and  in  the  boundary  conditions  for 
Eq.(2.8),  which  will  be  presented  shortly. 

Solution  of  Eq.(2.8)  with  appropriate  boundary  conditions 
provides  the  pressure  distribution  over  the  entire  region  o*  interest. 
Alternatively,  the  pressure  at  particular  points  in  a region  can  be 
obtained  by  integration  of  the  total  derivative  of  pressure  along  a 
given  path  [18].  If  the  pressure  is  known  at  point  A,  then  at  point  B 
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The  necessary  pressure  gradients  are  obtained  from  the  flow  field 
solution,  through  expressions  similar  to  Eqs.(2.1)  and  (2.2),  again 
illustrating  the  dependence  of  the  pressure  field  on  the  flow  field. 
This  approach  to  the  pressure  problem  is  popular,  for  example,  in 
lift  and  drag  calculations  for  a surface  where  only  the  pressure 
values  on  the  surface  itself  are  needed  [14,28]. 

Although  both  approaches  for  obtaining  pressure  values  would 
give  identical  results  in  the  continuum,  they  could  have  significant 
differences  in  the  discrete  case.  These  will  be  discussed  more 
fully  later  in  this  chapter. 

Solution  of  Eqs . (2.6) , (2.7),  and  (2.8)  will  model  the  flow 
through  any  geometry  when  subjected  to  appropriate  boundary  condi- 
tions. However,  with  simplification  of  the  boundary  conditions  and 
ease  of  programming  in  mind,  the  geometry  to  be  considered  was  modi- 
fied from  that  of  Fig.  1.1.  Modifications  include  the  following: 

a.  The  control  nozzle  and  control  ports  have  been  omitted. 
Their  effect  on  the  inlet  jet  is  modeled  as  a specified 
deflection,  cf> , of  the  jet.  The  addition  of  control  jet 
action  into  the  model  would  be  complex  with  regard  to  the 
programming  involved,  but  does  appear  to  be  theoretically 
straightforward. 
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b.  The  amplifier  considered  is  non-vented,  as  opposed  to  the 
vented  amplifier  of  Fig.  1.1.  This  is  realistic  in  the 
sense  that  non-vented  amplifiers  have  use  as  flow  divid- 
ers and  flow  signal  amplifiers.  However,  it  would 
greatly  expand  the  usefulness  of  the  model  if  a method 
for  handling  the  vent  flows  could  be  devised.  Such  a 
method  is  not  theoretically  available  at  the  present 
since  it  must  allow  for  unknown  flow  rates,  velocity 
profiles,  and  flow  reversals. 

c.  A square  cavity  region  is  considered.  This  is  solely  to 
simplify  the  programming  as  sloping  walls  could  be  han- 
dled if  desired. 

d.  The  center  vent  is  omitted.  Addition  of  the  center  vent 
would  be  possible  if  a flow  rate  could  be  somehow  estab- 
lished. A fully  developed  velocity  profile  corresponding 
to  this  flow  rate  could  then  be  assumed. 

e.  The  receivers  are  at  a 45°  angle  to  the  amplifier  center- 
line.  As  in  (c)  above,  this  is  for  programming  conven- 
ience only. 

The  resulting  geometry  is  shown  in  Fig.  2.1,  where  all 
lengths  are  non-dimensional  and  have  arbitrarily  selected  magnitudes. 
The  coordinate  system  is  fixed  such  that  the  origin  is  at  the  lower 
left  corner  of  the  cavity  region,  and  such  that  the  y-axis  is 
parallel  to  the  amplifier  centerline.  The  flow  enters  through  the 
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single  port  at  the  bottom  of  Fig.  2.1,  crosses  the  cavity  region, 
and  exits  through  the  two  receivers. 

B.  Boundary  Conditions 

For  the  geometry  of  Fig.  2.1,  all  boundaries  are  either  walls 
or  lines  on  which  something  is  known  about  the  velocity  profile  and 
flow  rate.  Continuum  boundary  conditions  for  these  cases  are  as 
follows. 

At  the  solid  walls  there  exists  a no  slip  condition.  Setting 
u and  v to  zero  in  Eq . (2.5)  the  no  slip  condition  implies  a constant 
ip  on  each  wall.  It  is  assumed  that  the  left  wall  of  the  amplifier 
has  i|>  = 0,  where  the  reference  value  0 was  arbitrarily  chosen.  As 
noted  earlier  the  flow  variables  are  non-dimensional i zed  with  respect 
to  the  average  inflow  velocity  and  the  supply  port  width.  Due  to 
this  choice  of  non-dimensional izing  variables,  the  right  wall  has 
^ = -1.0. 

The  remaining  solid  surface  is  that  positioned  between  the 
two  receivers,  which  will  be  referred  to  as  the  splitter.  The  value 
given  to  \p  on  the  splitter  will  be  used  to  model  different  loadings 
on  the  receivers.  As  in  all  flow  situations  the  flow  through  the 
receivers  is  proportional  to  the  pressure  drop  down  the  receiver. 

Thus,  changing  the  flow  by  changing  ip  on  the  splitter  will  model  a 
different  pressure  condition  at  the  outlet,  given  that  the  amplifier 
inlet  conditions  remain  the  same.  Then  for  a given  inlet  condition 
a characteristic  curve  of  receiver  flow  versus  pressure,  as  in 
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Fig.  1.5,  can  be  obtained  by  varying  the  flow  rate  in  the  receivers 
and  recording  the  resulting  pressures. 

The  new  outlet  pressure  is  not  known  until  after  the  flow 
solution  is  obtained.  Thus,  the  model  will  not  accept  a receiver 
outlet  pressure  as  a boundary  condition.  This  is  a result  of  the 
required  boundary  conditions  on  the  pressure  equation. 

At  the  inlet  a fully  developed  parabolic  velocity  profile  is 
assumed.  It  is  necessary,  however,  to  apply  this  condition  in  such 
a way  that  the  inlet  jet  can  have  an  initial  specified  deflection,  <!>. 
This  deflection  is  to  partially  model  the  effect  of  the  control  jets 
on  the  power  jet.  Although  the  control  jets  would  have  other  effects 
on  the  flow  field,  their  main  purpose  is  to  deflect  the  power  jet. 

As  pointed  out  earlier,  this  deflection  is  the  physical  mechanism 
which  allows  amplification  in  a fluidic  amplifier.  As  such  it  is 
important  to  incorporate  the  deflection  into  the  numerical  model. 

Consider  the  inlet  configuration  shown  in  Fig.  2.2.  With  n 
measured  perpendicular  to  left  wall  of  the  inlet,  a parabolic  velocity 
profile  leads  to  a stream  function  distribution  of 

3 2 

<Mn)  = 2 [ S)  ‘ 3 (w)  (2.10a) 

Since  it  is  desired  that  ip  vary  from  0 to  -1  across  the  inlet,  the 
value  of  w decreases  as  <|>  increases  in  magnitude,  but  the  inlet  width 
in  the  x direction  remains  constant  at  a value  of  1. 

Linder  a coordinate  transformation  from  the  (n,a)  to  the  (x,y) 
coordinate  system,  the  values  of  ip , c,  and  3tp/ 9x  remain  invariant. 
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Figure  2.2.  Amplifier  Inlet  Geometry 


Only  the  value  of  3ip/3y,  or  the  u velocity,  changes.  Since  3iJ;/3y 
was  the  only  variable  that  changed  with  deflection,  it  was  selected 
as  one  of  the  boundary  conditions  on  the  inlet  flow.  The  vorticity 
distribution  compatible  with  a parabolic  velocity  profile,  which  is 
derived  from  Eq.(2.10a),  was  selected  as  the  other  boundary  condition. 
These  two  boundary  conditions  are 

= -6  tan  <j>  [(x  - xQ)2  4 (x  - xQ) ] (2.10b) 

C = 12  (x  - xQ)  - 6 . (2.10c) 

At  the  receiver  outlets  the  flow  has  a developed  profile,  but 
the  restriction  of  a parabolic  profile  is  not  imposed  [20].  This 
assumption  yields  boundary  conditions  of  the  form 
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= 0 (2.11b) 

where  the  coordinate  axis  "a"  is  parallel  to  the  receiver  centerline. 

Switching  from  the  flow  field  boundary  conditions  to  consid- 
eration of  the  pressure  problem,  the  vorticity  and  stream  function 
equations  were  derived  in  such  a way  that  the  pressure  was  eliminated. 
It  follows  that  boundary  conditions  for  the  pressure  equation  must  be 
compatible  with  Eqs .(2.1)  through  (2.3).  The  result  is  gradient 
boundary  conditions  for. the  pressure  equation,  which  are  given  by 
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An  elliptic  partial  differential  equation,  such  as  Eq.(2.8), 
which  is  subject  to  gradient  boundary  conditions  on  all  boundaries, 
such  as  those  given  by  Eq.(2.12),  is  commonly  referred  to  as  a 
Neumann  problem.  For  this  type  of  problem  to  have  a solution  a 
further  constraint,  known  as  a consistency  equation,  must  be  met 
[8,15,29].  This  equation  gives  a unique  relationship  between  the 
gradients  of  the  dependent  variable  on  the  boundary  and  the  forcing 
function  of  the  partial  differential  equation.  For  the  pressure 
problem  this  relationship  is 
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where  n is  the  outward  directed  normal  from  the  boundary  and  s is 
tangential  to  the  boundary. 

A solution  to  Eq.(2.8)  which  meets  the  gradient  boundary  con- 
ditions given  by  Eq .(2.12)  is  not  a unique  solution  to  that  system  of 
equations  since  any  solution  plus  a constant  is  also  a solution  [8]. 

A solution  is  made  unique  when  it  is  also  desired  that  tnere  be  a 
fixed  value  for  pressure  at  some  point  within  the  boundaries.  In 
effect,  this  particular  value  can  be  used  to  fix  the  constant 
referred  to  above. 

C . Finite  Difference  Solution  for1  the  Flow  Field 

As  in  any  finite  difference  solution  to  a partial  differential 
equation,  the  general  procedure  here  is  to  discretize  the  vorticity 

transport  and  stream  function  equations  to  establish  a system  of 
algebraic  equations.  These  algebraic  equations  can  be  solved  for 
the  values  of  \p,  c,  u,  and  v at  all  points  on  the  finite  difference 
grid. 

The  grid  under  consideration  spans  the  region  bounded  by  the 

amplifier  boundaries  shown  in  Fig.  2.1.  The  region  is  divided  such 

that  in  the  x direction  there  are  I grid  points  separated  by  a 

m3  x 

distance  Ax.  Similarly,  in  the  y direction  there  are  J grid 

nlaX 

points  separated  by  a distance  Ay.  I w and  J were  chosen  to 

max  max 

produce  a grid  with  Ax  = Ay  = A,  although  this  is  not  necessary. 

The  notation  used  to  distinguish  a variable  value  at  a given  point 
is 
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0^  = 0 {(i  - 1)  Ax,  (j  - 1)  Ay}  (2.14) 

where  0 is  an  arbitrary  variable,  1 < i < I . and  1 < j < J,v. 

The  system  of  algebraic  equations  obtained  from  the  differ- 
ence equations  for  stream  function  and  vorticity  are  solved  using 
alternating  direction  implicit  (ADI)  methods.  The  methods  are  based 
on  the  work  of  Peaceman  and  Rachford  [17]  and  of  Douglas  [7].  Such 
methods  are  widely  used  in  finite  difference  work.  Descriptions  of 
the  method  can  be  found  in  many  books  on  finite  difference  methods, 
including  those  by  Roache  [20] , Mi tchel  1 [15],  Young  [29],  and  Young 
and  Gregory  [30].  Solutions  of  viscous  flow  problems  using  ADI 
methods  are  discussed  by  Roache  and  examples  are  found  in  papers  by 
Pearson  [18],  Briley  [6],  and  Richards  and  Chavez  [19].  A major 

reason  for  the  use  of  ADI  methods  is  that  the  splitting  of  directions 
produces  systems  of  equations  with  tridiagonal  coefficient  matrices, 
for  which  solutions  can  be  obtained  very  efficiently.  Another  reason 
for  its  use  is  its  excellent  stability  characteristics.  On  parabolic 
equations  this  stability  allows  ADI  methods  to  use  much  larger  com- 
putational time  steps  than  explicit  methods.  Thus,  ADI  methods  use 
less  computer  time  to  reach  a steady  state  solution.  Also,  the  con- 
vergence rate  of  the  method  is  good  in  solution  of  the  elliptic 
stream  function  equation  when  a sequence  of  iteration  parameters,  a 
p-sequence,  is  used. 

Along  with  selection  of  a solution  technique,  a differencing 
scheme  was  finally  selected  for  Eqs . (2.6)  and  (2.7).  Difficulties 
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with  stability  at  higher  Reynolds  numbers  led  to  the  use  of  various 
initial  differencing  schemes  on  the  nonlinear  terms  of  the  vorticity 
eguation.  In  particular,  second-order  accurate,  centered  differencing 
for  the  advection  terms  was  tried  and  abandoned  in  favor  of  the  more 
stable,  first-order  accurate,  first  upwind  differencing  [20].  This 
method  proved  to  have  convergence  problems  at  regions  where  velocity 
reversals  occurred  and  was  replaced  by  second  upwind  differencing 
[20].  Although  different  in  algebraic  form,  this  difference  scheme 
can  be  shown  to  be  the  same  as  that  used  by  Runchal , Wolfshtein  and 
Spalding  [23]. 

Centered  differencing  is  used  on  the  diffusion  terms.  The 
resulting  difference  equations  for  use  with  ADI  in  solving  the  vor- 
ticity transport  equation  are 
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With  the  linear  stream  function  equation  second-order  accu- 
rate, centered  differencing  was  used  on  both  second  derivatives. 

The  resulting  difference  equations  for  use  with  ADI  are 
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x-sweep: 


Cj  - <2  + ».>  AT  * AT,S 

- - - (2  + Pt)  »;j  + - *x2c(j 


(2.16a) 


y- sweep: 


<]♦,  • (2  + “*>  + <]-i 

■ - [*}3j  - <2  * p.)  C * ATj  - 


Ay2?ij 


(2.16b) 


where 


k = iteration  index 
k+^  = intermediate  computation  step 
= iteration  parameter,  1 < H < m . 

The  optimum  sequence  of  iteration  parameters  p„  is  not  known 
for  the  problem  being  considered.  However,  a sequence  based  on  the 
Wachspress  parameters  [29]  worked  quite  well.  For  purposes  of 
evaluating  this  sequence  the  region  of  interest  was  taken  as  the 
largest  square  in  the  complete  amplifier.  This  corresponds  to  the 
cavity  region  and  is  a square  spanned  by  a grid  of  Imax  by  Imgx 
points.  The  Wachspress  parameters  are  given  by 

ft-1 

p^  = 3 (a/3)m_1  , m > 2,  SL  = 1,  2,  ...,  m (2.16c) 
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where  a and  B can  be  estimated  by 


. . 2 IT 

a = 4 sin 


2 I 


max 


B = 4 cos^ 


(2. 16d) 

(2.1 6e ) 


max 


and  m is  the  smallest  integer  such  that 


(0.414) 


2m 


< a/ 6 


(2. 16f) 


For  I. 


max 


41 , this  resulted  in 


and 


m = 4 


p£  = 3.9939  (0.001469) 


fc-1 

3 


Z = 1 , 2,  3,  4 


(2. 16g) 


Equations  (2.15)  and  (2.16)  will  pro1  : ie  the  necessary 
algebraic  equations  at  all  points  in  the  interior  of  the  grid.  How- 
ever, at  the  boundaries  they  would  require  information  from  points 
outside  of  the  grid.  Here  the  information  contained  in  the  continuum 
boundary  conditions  is  used  to  form  the  necessary  discrete  equations. 

On  the  solid  walls  the  discrete  values  of  ip  are  simply  the 
continuum  values  of  <Jj.  Then,  for  nodes  on  the  walls 


i »J 


0.0, 


spl itter’ 


left  wall 
splitter 


(2.17) 


-1 .0,  right  wal 1 

Two  separate  equations  are  used  for  the  boundary  values  of 
vorticity  on  the  walls,  one  in  the  cavity  region  and  one  in  the 
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receivers.  Both  are  derived  using  the  no  slip  wall  condition.  For 
the  cavity  region  the  derivation  involves  approximating  with  a 
four  point,  cubic,  LaGrangian  interpolating  polynomial.  This  can  be 
differentiated  twice  arid,  using  the  velocity  on  the  wall  and  tne 
definition  of  vorticity,  evaluated  at  the  wall.  The  result  is  a 
second-order  accurate  approximation  for  r,  on  the  wall  [6], 

For  example,  consider  the  wall  point  shown  in  Fig.  2.3. 


® (1,3) 

<g>  (i,2) 

-7-7" X -7  / / 7 /S>/[r'7~Jy 7 * / / / / 

® (1,0) 

Figure  2.3.  Grid  Notation 

For  these  points 

^(y)  = ^3  [-  (y-yl)(y-y2)(y-y3}  ^i.O 

+ 3 (y-y0)(y-y2)(y-y3) 

- 3 (y -yQ)(y-y1)(y-y3)  ^i>2 
+ (y - y0)(y - )(y - y2)  ^ 


(2.18a) 


1 
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where  y^  denotes  y at  the  point  (i,j).  Differentiating  this  once  and 
evaluating  on  the  wall  gives 


U- 


i ,1 


3y 


= 0 = 


1,1 


= Ay  ('  2^i,0  ‘ 3*i,l 


+ 6l^i,2  ' ^i,3)  • 


Differentiating  again  and  evaluating  at  the  wall  gives 

-V  <*1i0  - 2*i,,  * *,,2>  ■ 


2. 
3 ip 

3? 


, ,1 


Ay 


Since  <p  is  constant  on  the  wall 


_ 

^ ' V 


i ,1 


(2.18b) 


(2.18c) 


(2.1 8d ) 


Eliminating  ip.  n from  Eqs. (2.18b)  and  (2.18c)  and  substituting  with 
1 , u 

Eg. (2. 18d) , 


51.1  3 7-7  (-  7^i,l  + 8*i,2  ' ■ 


(2.1 8e ) 


2 Ay 


As  recommended  by  Briley  [6],  to  increase  stability  the 
velocity  at  the  point  adjacent  to  the  wall  is  also  evaluated  from 
Eq. (2.18a).  The  point  (i,2)  in  Fig.  2.3  is  such  a Doint.  Its 
u-velocity  is  given  by 


u • 


i,2 


4Ay 


(-  5iPU}  + 40i>2  + ^ij3) 


(2. 18f ) 


The  values  of  a on  the  sloping  walls  of  the  receivers  can  be 
derived  using  two  Taylor  series  expansions  for  'p  at  a wall  node.  For 


' 
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example,  consider  the  left  wall  of  the  left  receiver.  About  a point 
(i  ,j)  on  the  wall 

,2., 


j s *1  ,j f S 


*.  * *.  . + |s 

Vj+l  vi,j  Jy 


Ax  + bAx2  ^ 
i,j  3x 


Ay  + ^Ay2  Mr 

i.j  9y 


+ 0(Ax°)  (2.19a) 


i.J 


+ 0(AyJ).  (2.19b) 


i.J 


Due  to  the  no  slip  condition  the  first  derivatives  are  zero  and  can 
be  eliminated.  Further  simplification  is  achieved  by  replacing  Ax 
and  Ay  with  A and  adding  the  two  equations.  Rearranging  and  using 
Eq.(2.7),  the  expression  for  vorticity  is 


2 


Ci,j  ' 7(*vH,j  +*1,j+l  ' + °(6)  • 


(2.19c) 


Similar  expressions  hold  for  t;.  . on  the  other  receiver  walls. 

■ *J 

In  the  inlet  the  vorticity  boundary  condition  is  given  by  the 
discrete  analog  of  Eq. (2.10c),  which  is 


’i  .1 


12  ((i  - 1)  Ax  - xn)  - 6 


(2.20) 


The  boundary  condition  on  ^ in  the  inlet  is  derived  from  a Taylor 
series  expansion  for  ip  about  a point  in  the  inlet,  that  is 

+ Ay^  3^  +0(Ay3)  . (2.21a) 

1 ^ i,l 

The  first  derivative  in  this  expression  is  obtained  from  Eq.(2.10b). 

An  expression  for  the  second  derivative  is  obtained  from  Eq.(2.7), 

2 2 

where  <;•  i is  known  and  3 iJj/3x  can  be  evaluated  using  centered 
1 * ' 


• *i,i  + 
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differences.  Substituting  into  Eq. (2.21a)  the  resulting  expression 
is 

p 

ip..  i = h (2iJk  2 - 12Ay  tan  <j>  { ( i Ax  - Ax  - xQ) 

+ iAx  - Ax  - xQ}  - ^2Cij  + ^i+i  j (2.21b) 

+ Vt  ,i^  + °(Ay2) 


which  would  be  second-order  accurate  due  to  the  differencing  of 

? ? 

3 ip/3x  . 

At  the  receiver  outlets  one-sided  difference  expressions 
were  used  for  the  first  derivatives  in  Eq .(2.11).  For  the  two  re- 
ceivers this  resulted  in  four  expressions: 

Left  receiver, 


?i,j  = Ci+l,j-l 

, j = ^i+l.j-l 


(2.22) 


Right  receiver, 

Ji,j  " ?i-l,j-l 


i-1 ,j-l 


(2.22) 


Special  care  was  given  to  the  values  of  vorticity  at  the  sharp 
corners  in  the  amplifier  geometry.  At  those  six  corners  which  are 
concave  the-'e  is  no  problem,  since  these  points  never  enter  into  the 
calculations.  However,  at  the  convex  corners  of  the  receiver  inlets 
the  proper  value  of  c to  use  was  in  question.  Following  Roache's 
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recommendations  [20],  the  procedure  being  used  is  to  allow  r to  be 
discontinuous  at  these  points.  If  the  value  of  t,  is  needed  for  cal- 
culations dealing  with  mesh  points  inside  the  receivers,  r is  defined 
as  in  Eq. (2.19c).  If  the  flow  outside  the  receiver  in  the  main 
amplifier  cavity  region  is  being  considered,  £ is  defined  as  in 
Eq. ( 2 . 1 Be ) . 

At  the  receiver  outlets  one  corner  point  enters  the  calcula- 
tions. However,  its  value  is  defined  by  considering  the  point  to 
have  an  extension  of  the  wall  beyond  it.  Thus  c at  this  point  is 
evaluated  as  in  Eq . ( 2 . 22 ) . 

At  the  corners  of  the  amplifier  inlet  the  value  of  vorticity 
is  again  discontinuous.  The  two  choices  considered  for  defining  this 
value  were  either  to  base  it  on  the  values  of  <ii  in  the  adjacent  flow 
field,  as  in  Eq. (2.19c),  or  to  include  it  as  part  of  the  inlet  pro- 
file, as  in  Eq. (2. 20) . Based  on  the  comparison  of  numerical  and  ex- 
perimental results,  vorticity  for  these  two  points  was  defined  as  in 
Eq. (2.20) . 

Using  an  implicit  method  on  the  vorticity  equation  can  pro- 
duce a problem  with  the  boundary  conditions.  The  solution  for  the 
n+1  time  step  requires  boundary  values  which  are  not  available  [20], 
since  they  also  need  to  be  at  the  n+1  level.  Further,  using  lagged 

boundary  conditions,  i.e.  using  boundary  conditions  based  on  informa- 
t h 

tion  at  the  ncn  time  level,  has  been  shown  to  be  numerically  desta- 
bilizing [20],  The  procedure  chosen  to  resolve  this  problem  is  to: 
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a.  Extrapolate  using  past  values  of  r,  for  an  estimate  of 

w 

rn  + \ Here  the  subscript  w denotes  a wall  value. 

'w 

b.  Extrapolate  using  past  values  of  u,  ■ and  v.  . for  an 

• » J ' > J 

estimate  of  u?+]  and  v^+.. 

c.  Calculate  new  values  of 

• »J 

d.  Calculate  new  values  of 

e.  Calculate  new  estimates  for  uV+] , v?+]  and  tn+1  using 
. n+1 

f.  Compute  the  change  in  wall  vorticity. 

g.  If  convergence  occurs,  as  will  be  discussed  later,  go 
to  the  next  time  step.  If  not,  go  back  through  (c)  to 
(g)  until  convergence  is  achieved  or  until  a prescribed 
number  of  iterations  have  been  completed. 

This  procedure  is  suggested  by  Roache  and  is  similar  to  that  used  by 
Briley  [6]. 

It  is  also  possible  with  the  above  iteration  scheme  to  set  a 
variable  time  step  based  on  the  numbers  of  iterations  for  convergence 

of  the  previous  time  step  [6],  This  allows  the  solution  to  use  a 
time  step  which  is,  in  a sense,  numerically  optimal. 

There  are  three  convergence  criteria  involved  in  this  solu- 
tion technique.  These  are  convergence  of  the  overall  flow  field, 
convergence  of  the  boundary  values  of  vorticity  at  each  time  step, 
and  convergence  of  the  stream  function  at  each  time  step.  The  over- 
all solution  is  considered  converged  when  the  normalized,  relative 
change  in  vorticity  over  a time  step  is  less  than  some  c for  each 
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point  in  the  mesh.  That  is,  the  solution  is  converged  when  for  all 
points 


i n+1  .n 

, kiax  knax 

1 r r [ i 

1 

At 

k - a < 

np  ■-  £ ki  i 

i=l  j=l  1 1,J|j 

X 

] 

At  , 

L ref 

where 

NP  = number  of  points  in  the  flow  field 

At  c = reference  time  step, 
ref  r 

A typical  value  for  c is  0.001. 

The  convergence  criteria  on  the  boundary  values  of  vorticity 

and  on  the  stream  function  are  interrelated  and  designed  so  that  they 

become  more  stringent  as  the  overall  solution  converges.  As  the  flow 

field  solution  approaches  steady  state  the  change  in  vorticity  across 

a time  step,  cn+1  - tn|,  becomes  smaller.  Therefore,  the  allowed 

W to 

change  between  estimates  of  must  be  reduced.  To  implement  this 
check  on  convergence,  the  criterion  used  is  that  at  each  wall  point 

I » n ^ I I ..  . n _ I 


LCP  L 


kn+1-cn 
w -w 


WZTEST 


(2.24) 


where 


q = an  iteration  counter 
LCP  = number  of  boundary  points 

= constant  factor,  typically  0.5. 

This  criterion  prevents  unnecessary  computational  labor  at  early 
solution  times  but  guarantees  that  convergence  will  always  be  suffi- 
ciently stringent. 
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Also,  since  wall  vorticity  is  calculated  from  stream  function 
values  at  adjacent  grid  points,  Eqs. (2.19c)  and  (2.18e),  the  stream 
function  convergence  must  be  related  to  the  criterion  for  wall  vor- 
ticity convergence.  Upon  acceptable  convergence  of  the  stream  func- 

I k+1  k 1 

tion  the  allowed  magnitude  of  change,  (K  • - • , must  be  less 

| ' ' »J  I 

than  the  magnitude  of  change,  WZTEST,  allowed  for  the  wall  vorticity. 

As  an  example,  consider  the  following.  If  the  convergence  on  \p  were 
k + 1 k 

such  that  iji.  . - tjj.  . = 0.001,  then  the  variation  in  r in  Eq.  (2.19c) 

i .]  i J w 

due  solely  to  lack  of  convergence  of  ip  would  be  on  the  order  of 

p 

0.002/Ax  . This  could  easily  be  greater  than  the  allowed  change  in 
Cw-  To  prevent  such  an  occurrence  the  stream  function  at  each  grid 
point  is  converged  to 

< F,  x WZTEST  (2.25) 

i J - 2 

where  F0  is  a constant  factor,  typically  0.005. 

This  completes  the  discussion  of  the  discrete  solution  tech- 
nique for  the  flow  field.  In  summary  the  solution  sequence  is: 

a.  Assume  an  initial  distribution  for  • , the 

necessary  boundary  conditions,  and  velocities. 

b.  Calculate  using  ADI. 

c.  Calculate  i|in+1  using  ADI. 

d.  Calculate  new  boundary  conditions  for  ?n+  . 

e.  Iterate  on  (b)  through  (d)  until  convergence  criterion 

on  r"  is  satisfied, 
w 
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f. 


g- 


h. 


i . 


Check  overall  convergence  of  -.  Stop  if  solution  con- 
verged. 

Update  the  time  step  index,  varying  the  time  step  if 
desired. 

Extrapolate  for  and  vr.l+|  velocity  components  and 

n+1  , 

r values, 
w 

Return  to  (b). 


D . Method  of  Solution  for  the  Pressure  Field 

As  pointed  out  earlier,  two  separate  approaches  could  be  taken 
in  solving  for  the  necessary  pressure  values.  The  method  involving 
integration  along  a path,  Eq.(2.9),  has  at  least  one  possibly  serious 
disadvantage.  Since  the  static  pressure  is  a point  function  and  not 
a path  function,  integration  between  points  A and  B on  two  separate 
paths  should  produce  the  same  value  for  the  pressure  at  B.  However, 
as  Roache  [20]  points  out  this  could  be  very  difficult  to  achieve 

numerically.  In  obtaining  the  pressure  gradients  to  be  integrated, 
error  would  be  introduced  by  the  numerical  differentiation  necessary 
to  evaluate  Eq . (2.12).  Error  would  also  be  introduced  by  the  numeri- 
cal integration  of  the  resulting  values. 

Because  of  these  problems,  Roache  recommends  the  solution  of 
Eq.(2.8)  for  the  desired  pressure  field.  That  approach  was  adopted, 
partially  on  this  recommendation  and  Dartially  because  solution  of 
Eq.(2.8)  was  anticipated  to  be  straightforward.  Except  for  the 
change  in  boundary  conditions  this  is  the  same  type  of  system  which 
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governs  the  stream  function,  Eq .(2.7).  However,  solution  of  the 
pressure  system  has  proven  to  be  considerably  more  difficult,  due 
both  to  the  gradient  boundary  conditions  and  to  the  evaluation  of 
the  forcing  function.  In  particular,  the  calculations  are  compli- 
cated by  a failure  of  the  system  to  satisfy  the  discrete  analog  of 
Eq . (2.13) , the  consistency  equation. 

As  for  the  flow  field,  a set  of  simultaneous  algebraic  equa- 
tions must  be  established  from  the  discrete  analog  to  Eq.(2.8).  This 
set  is  of  the  form 

A P = S 

where  if 

N is  the  number  of  unknown  pressure  values, 

then 

P is  a (Nxl)  vector  of  the  pressures  P.^, 

A is  a (NxN)  matrix  of  difference  coefficients,  arid 
is  a (Nxl)  vector  containing  values  for  the  right-hand 
side  of  Eq.(2.8)  and  for  the  gradient  boundary  conditions. 

As  a check  on  the  differencing  techniques  used,  Eq.(2.26) 
must  reflect  the  fact  that  the  Neumann  problem  of  Eqs.(2.8)  and 
(2.12)  has  no  unique  solution.  However,  because  a unique  solution 
may  be  found  in  terms  of  one  given  pressure  value  at  one  point  in  the 
field,  it  has  to  follow  that  A in  Eq.(2.26)  is  singular  and  of  rank 
N-l . 

Expanding  this  idea,  consider  each  column  of  A,  denoted  as 

a , as  a vector  in  an  N-space.  Define  an  additional  vector  c , 

-m  - 
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cT  - [cj , c2,  c3,  cm,  cN]  • (2.27) 

If  a c can  be  found  whicfi  is  mutually  orthogonal  to  each  a^,  then 

there  has  to  be  at  least  one  linear  dependency  in  the  set  of  vectors 

. (Only  N vectors  in  an  N-space  can  be  mutually  orthogonal.) 

Further,  such  a vector  will  also  be  orthogonal  to  the  S vector.  If 

it  is  assumed  that  there  is  only  one  linear  dependence  in  the  set 

a , it  can  be  concluded  that  A is  of  rank  (N-l).  Thus,  the  differ- 
-m 

encing  technique  used  would  retain  the  characteristics  of  the  Neumann 
problem. 

Practically,  for  a proper  difference  formulation  of  the 
Neumann  problem,  a matrix  C can  be  found  by  inspection  such  that 

C(NxN)  = dia9  (C1’  C2  ’ c3 V CN}  (2'28) 

Note  that  due  to  the  definition  of  c_  the  columns  of  the  product  CA 
will  sum  to  zero  and  the  elements  of  the  vector  CS  will  sum  to  zero. 
This  second  summation  is,  in  fact,  equivalent  to  a numerical  evalua- 
tion of  Eq . (2.13)  with  the  proper  choice  of  integration  formulas. 

This  summation  will  be  referred  to  as  the  discrete  analog  of  the  con- 
sistency equation.  The  above  information  provides  a method  for  eval- 
uating and  developing  the  necessary  difference  equations. 

In  particular,  the  summation  of  elements  in  the  columns  of 
CA  is  a check  on  the  implementation  of  boundary  conditions.  If  not 
properly  implemented,  the  summations  cannot  all  be  made  equal  to  zero. 

With  this  in  mind,  difference  approximations  were  derived  for 
the  interior  and  boundary  grid  points.  For  interior  grid  points  the 
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standard  five  point  centered  difference  approximation  was  used  to 
difference  Eq.(2.8),  that  is 


P. 


i ,J 


^ ( P . . -f  P . _ . + P . . _ + P . . 

i+l  ,J  i-l ,J  i ,j  + l i ,J-1 


(2.29) 


where 


ro 

,2.  , 

,2.  , 

2 | 

6 Ip 

— 7 

o ip 

— 2 ! 

• • 6y  • . 

i.j  h.j 

6 t >i 

li.jj 

J 

At  non-corner  points  of  the  cavity  boundary,  the  discrete 
boundary  equation  is  similar  to  Eq. (2.21b)  and  is  derived  in  the 
same  manner.  For  instance,  on  the  lower  cavity  boundary 


P 


i ,1 


P. 


i + l ,1 


+ P 


i-l  ,1 


...  <SP 
Wlth  6y 


I i >1 


obtained  by  differencing  Eq . (2.12). 


(2.30) 


The  difference  equation  is  similar  for  the  four  concave  cor- 
ner points  on  the  boundary  in  the  cavity  region  and  for  the  non- 
corner points  on  the  receiver  walls.  As  shown  in  Fig.  2.4,  in  both 
cases  only  two  adjoining  grid  points  are  available  to  form  the  needed 
difference  approximation. 


r 
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a)  Lower  Left 
Corner 

b)  Upper  Left  c) 

Corner 

Upper  Right  d) 
Corner 

Lower  Right 
Corner 

<8> 

<8) 

<8> 

/-ri^7-7-7&kr-r 

e ) Left  Wall, 

Left  Receiver 

\ 

f ) Right  Wall , g) 
Left  Receiver 

Left  Wall , h) 

Right  Receiver 

^ <g> 

Right  Wall, 
Right  Receiver 
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i)  Corner,  Left  j) 
Receiver 

Corner,  Right 
Receiver 

AT  -&■ - 0— 

A 

A 

1 

<S> 

1 

\ 

/ 

Figure  2.4.  Grid  Configurations  for  Concave 
Corners  and  Receiver  Walls 


Let  the  point  being  considered  be  (i,j).  Let  the  adjoining  point  in 
the  x-direction  have  coordinates  (i',j),  where  ' = i 1 . Similarly, 
let  the  adjoining  point  in  the  y-direction  be  at  ( i , j 1 ) , where 


j'  = j±l.  Define 

ID  as 

(i'-i) 

Ax  and  JO 

as 

(j’  - J ) Ay . 

Then,  using 

Taylor  series  expansions 

pi'j  ■ 

P.  . + 
i ,J 

ID  !x- 

AX2 

,.i+^ 

a2o 

9X2 

+ 0(.ix3) 

(2.31a) 

pi,i'  * 

P.  . + 
i ,J 

JD  fy 

.*4 

I'.O 

. 2 
3y 

+ 0(Ay3)  . 

(2.31b) 
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ing  and  using  Eq . (2-8)  to  introduce  s.  ., 

' >J 


P.  . = H 


P.,  . + P. 
i ,J  i,j 


ID  ~ 

dX 


JD 


i ,J 


9P 

ay 


(2.31c) 


~2~  Si ,jj 


+ 0(AJ)  . 


At  the  inlets  and  outlets  to  the  receivers  the  difference 
formulas  had  to  be  matched  with  the  surrounding  difference  patterns 
in  order  to  insure  the  singularity  of  A.  As  pointed  out,  the  criter- 
ion for  this  was  the  summation  of  the  columns  of  CA.  The  equations 
were  derived  by  considering  the  derivatives  of  pressure  discontin- 
uous, writing  Taylor  series,  and  summing  these  series  in  such  a way 
that  the  criterion  was  met.  The  derivation  for  the  point  of  the 
splitter,  point  2 on  Fig.  2.1,  will  be  presented.  The  other  deriva- 
tions are  similar  and  only  the  results  will  be  presented.  At  the 
splitter  point 


P, 


i-1  ,j 


Pi+1  ,J 


P.  . 
i ,J 


P.  . + Ax 


dp 

9x 


+ 0(Ax2) 
+ 0(Ay2) 
+ 0(Ax2) 


(2.32a) 

(2.32b) 

(2.32c) 


The  subscripts  on  the  derivatives  denote  the  coordinate  direction  in 
which  the  discrete  derivative  is  evaluated.  If  then  variables 
at  points  to  the  left  or  points  below  (i,j)  are  used  in  the  deriva- 
tive evaluation.  If  "+",  the  variables  to  the  right  or  above  (i,j) 
are  used.  Adding  the  above  equations 
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For  points  1 and  3 on  Fig.  2.1,  the  results  are 


For  1: 
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For  3: 
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(2.34) 


As  illustrated  in  Fig.  2.4,  the  grid  points  on  the  acute 
angle  corners  of  the  receiver  outlets  are  omitted.  Both  corners  of 
each  receiver  are  then  differenced  as  in  Eq. (2.29c).  The  remaining 
outlet  points  were  differenced  as 


P =’A'P  +P  +P  + a IE 

i.j  h 1/iJ-l  Ki+lJ  i-l,j  ay 


+ O(Ay^) 


(2.35) 
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The  system  of  algebraic  equations  formed  from  the  above 
difference  equations  could  be  solved  with  either  ADI  or  point  succes- 
sive over-relaxation,  SOR  [20,29].  SOR  was  chosen  because  of  its 
simplicity  and  the  ease  with  which  it  can  be  applied  and  because  the 
speed  of  convergence  was  not  of  prime  importance.  For  this  problem 
the  SOR  algorithm  can  be  written  as 


Pk+!  = (1  - uyj  Pk  . + wnP*  . 

i ,3  0'  i ,.]  0 i ,j 


(2.36) 


Here  u)„  is  the  relaxation  factor,  k is  the  iteration  index,  and  P.  . 

U i > J 

is  the  value  of  P.  . given  by  the  appropriate  equation  from  the 
group  of  Eqs.(2.27)  through  (2.33). 

In  general  there  exists  an  optimum  value  of  the  factor  u>, 

I < < 2.  However,  for  the  Neumann  problem  on  a non-rectangular 

grid  there  is  no  analytic  expression  for  this  value.  To  estimate 
the  procedure  adopted  was  to  calculate  Wq  for  use  with  Poisson's 
equation  and  Dirchilet  boundary  conditions.  As  discussed  by  Roache 
[20]  a value  slightly  larger  than  this  was  then  used  for  Wq.  For  a 
41  x 57  grid  the  u>q  chosen  was  1.9 

The  evaluation  of  s.  . , (6P/6y)|.  • and  (6P/6x)|.  • remains 
to  be  discussed.  These  values  must  be  obtained  through  numerical 
differentiation  of  the  discrete  flow  field  solutions.  The  complete 
set  of  difference  equations  is  presented  in  Table  A.l  of  the  Appendix. 
There  are,  however,  several  general  comments  that  need  to  be  made. 

The  numerical  differentiation  is  complicated  by  the  need  for 
one  sided  differencing  on  solid  walls  and  away  from  discontinuities 


ja. 
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in  vorticity.  This  problem  is  helped  somewhat  by  narrowing  the 
region  of  interest.  The  pressure  field  of  interest  is  that  near  the 
jet  as  it  issues  from  the  inlet,  flows  across  the  cavity,  and  flows 
down  the  receivers.  Since  the  pressure  gradients,  and  thus  the 
boundary  conditions,  can  be  evaluated  at  any  point  and  since  there 
is  no  need  for  the  pressure  values  in  the  cavity  side  regions,  only 
the  central  half  of  the  cavity  region  is  included  in  the  pressure 
solution.  This  results  in  the  new  boundaries  shown  by  the  dotted 
lines  in  Fig.  2.1.  Besides  reducing  the  total  number  of  SOR  calcula- 
tions, this  allows  the  use  of  central  differencing  for  both  s.  . and 

* 

( SP/c)y ) j . . on  the  new  boundaries. 

The  second  comment  concerns  the  evaluation  of  derivatives  in 
the  receivers.  In  this  region  of  relatively  high  gradients  normal 
to  the  wall  there  are  relatively  few  grid  points.  Thus,  in  the  nor- 
mal discrete  x-y  grid  system  the  values  of  s,  . are  noticeably  in- 

' >J 

accurate.  For  example,  towards  the  end  of  the  receivers  the  flow 
becomes  fully  developed.  For  this  case,  s.  ^ should  be  zero.  This 
result  was  not  obtained  with  differencing  in  the  normal  coordinate 
system.  For  this  reason  a transformation  of  coordinates  was  made  in 
each  receiver  such  that  the  new  system  had  one  axis  parallel  to  the 
centerline  of  the  receiver.  Such  a coordinate  system  is  shown  in  the 
left  receiver  of  Fig.  2.1.  This  transformation  introduced  discrete 
derivatives  in  the  axial  direction.  These  derivatives  go  to  zero 
as  is  required.  For  instance, 

f 


Since  the  streaml  ines  become  parai  lei  to  the  walls  at  the  end  of  the 
receivers,  the  ^ values  do  not  change  in  the  axial  direction.  There- 
fore, 6vb/6a  and  go  to  zero  and  s,  • goes  to  zero. 

The  third  comment,  and  certainly  the  most  important,  is  that 
no  difference  scheme  was  found  that  would  produce  values  for  the 
gradients  and  the  forcing  functions,  s — , that  would  satisfy  the 

discrete  consistency  equation.  A discussion  of  what  this  means  physi- 
cally, what  the  effects  on  the  numerical  solution  are,  and  possible 
causes  will  be  reserved  for  later.  At  this  point  a method  for  compen- 
sating for  the  inconsistency  difficulty  will  be  presented. 

It  was  noted  during  the  study  that  for  inconsistent  formula- 
tions, i.e.  those  for  which  the  elements  ot  CS  did  not  sum  to  zero, 
the  SOR  algorithm  would  converge  to  values  for  which  v2P|  . . / s-  .. 
For  every  point  the  solution  was  in  error  by  the  same  amount,  s', 
where  s‘  = V PI-  s,  _ . Further,  the  value  of  the  numerical 

| * J J I jJ 

integral  oi  o n , the  re<  i as  equal  to  the  summation  of  the 
elements  of  CS,  the  is,  to  the  inconsistency.  Therefore,  weighting 
the  contribution  of  each  point  by  its  contribution  to  the  numerical 
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integral  and  knowing  the  total  inconsistency,  s'  is  a calculable 
quantity. 

The  method  for  compensating  for  inconsistency  is  based  on  this 
calculation  of  s'.  After  evaluation  of  the  boundary  gradients  and 
s.j  . the  inconsistency  is  calculated.  From  the  inconsistency,  s'  is 
determined.  This  value  is  then  combined  with  each  s.  .in  such  a 
way,  meaning  either  added  or  subtracted,  that  the  inconsistency  is 
forced  to  zero. 

This  compensation  is  not  proposed  as  a correction  for  incon- 
sistency or  as  a completely  satisfactory  solution  to  the  problem.  It 
is  simply  a means  of  eliminating  some  of  the  obviously  incorrect  re- 
sults produced  by  inconsistency.  There  are  several  noticeable  effects 
of  inconsistency  which  will  be  discussed  in  detail  later.  Briefly  an 
inconsistent  formulation  results  in  solutions  which 

a.  Do  not  have  the  gradients  specified  by  the  boundary 
conditions. 

b.  Do  not.  produce  the  specified  forcing  function. 

c.  Do  not  converge  to  symmetric  pressure  fields  for 
symmetric  flow  fields. 

Although  not  necessarily  producing  more  accurate  solutions,  the  incon- 


sistency compensation  does  alleviate  these  problems. 

Solution  of  the  pressure  equation  subject  to  all  gradient 


I 
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this  shape,  not  to  some  unique  set  of  pressure  values.  The  following 
procedure  was  used  to  check  for  this  characteristic: 

a.  Select  one  point  in  the  grid  as  a reference  point,  say 
(g,h),  and  select  some  reference  value  of  pressure  for 
this  point,  say  1 .0. 


b. 

c. 

d. 


Determine  the  deviation  of  P . from  its  reference  value. 

9 >n 

Translate  the  entire  surface  by  adding  the  deviation  to 

each  P..  y Pg  ^ then  equals  the  reference  value. 

Calculate  new  values  of  P-  • using  SOR  and  allowing  P , 

i * j y »** 

to  change.  Then  translate  the  new  pressure  surface. 

nk  + l nk  n 

Now  P„  , - P„  h = 0. 
g»h  g,h 

k+1  k 

Determine  the  largest  deviation  between  P.  . and  P.  ., 

i.J  i.J 


i .e. 


pk+1 

i .j 


pjf  . 

i ,j | max 


As  the  surface  shape  converges 


this  deviation  will  approach  zero, 
f.  Terminate  the  calculations  when 


pk+1  _ pk 

i » J i ,j 


max 


(2.38) 


where  e is  a fixed  number,  typically  10  . 

The  above  procedure  is  similar  to  that  proposed  by  Fanning  and 
Mueller  [8], 

This  completes  the  discussion  of  the  solution  technique  for 
the  pressure  field.  In  summary  the  solution  sequence  is: 

a.  Calculate  s.  the  necessary  gradient  boundary  conditions, 

' ij 

and  compensate  for  inconsistency. 

b.  Assume  an  initial  distribution  for  pressure. 
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J( 

Translate  the  P.  ■ values. 

• i J |(4  ] 

Calculate  Pk+1  using  SOR.  Note  that  new  values  of  Pi j 

on  the  boundary  are  also  calculated. 
k+1 

Translate  the  P_  • values. 

Check  for  convergence  of  the  surface.  Stop  if  solution 
converged. 

Store  Pk+1.  in  Pk  • array  and  return  to  (d). 
i 1 >3 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

This  chapter  discusses  the  objective,  equipment,  and  proce- 
dures for  the  experimental  portion  of  this  study.  This  includes  a 
detailed  description  of  the  design  and  construction  of  the  experi- 
mental model  itself  and  of  the  associated  fluid  supply  system.  Also 
included  in  this  chapter  is  a discussion  of  the  instrumentation,  dye 
injection  techniques,  and  the  photographic  techniques  which  were 
used  for  the  flow  visualization. 

A.  Overview 

The  purpose  of  the  experimentation  was  to  provide  data  for 
verifying  the  numerical  simulations.  Basically,  a model  was  designed 
and  constructed  such  that  the  flow  in  the  horizontal  midplane  is  a 
laminar,  essentially  two  dimensional  flow  and  is  bounded  by  the  geom- 
etry shown  in  Fig.  2.1.  Streamline  patterns  were  produced  by  inject- 
ing dye  into  this  midplane.  These  patterns  can  be  qualitatively 
compared  to  the  numerically  predicted  streamline  patterns  for  the 
same  flow  conditions.  This  flow  visualization  is  used  as  qualita- 
tive verification  of  the  predicated  flow  fields. 

Because  of  practical  complexity  in  constructing  a model  to 
allow  deflection  of  the  power  jet,  only  the  symmetric  undeflected 
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case  was  studied  experimentally.  It  is  felt  that  pattern  verification 
in  this  case,  which  can  include  asymmetric  loading,  would  be  suffi- 
cient. 

Specific  design  requirements  of  the  model  and  experimental 
apparatus  were  as  follows: 

a.  The  geometry  of  the  test  section  must  be  identical  to 
that  of  Fig.  2.1.  Model  walls  must  be  transparent  to 
permit  flow  visualization. 

b.  The  flow  in  the  test  section  must  be  laminar  and  free  of 
small  amplitude  disturbances. 

c.  Inlet  and  outlet  flow  rates  must  be  adjustable  and  must 
be  measured. 

d.  Dye  must  be  injected  so  as  to  illustrate  the  flow  pat- 
terns but  not  disturb  the  flow. 

e.  Photographic  equipment  and  lighting  must  be  provided  for 
recording  the  flow  patterns. 

These  design  requirements  led  to  the  experimental  apparatus 
shown  in  Fig.  3.1.  This  apparatus  can  be  divided  into  the  model,  the 
fluid  supply  system,  instrumentation,  and  photographic  equipment. 

B.  Model 

The  model  designed  to  meet  these  requirements  is  fairly  sim- 
ple in  design  but  was  rather  difficult  to  construct.  A plan  view  is 
shown  in  Fig.  3.2.  The  walls,  cover  plate  and  base  plate  were  made 
of  clear  plexiglass.  Use  of  plexiglass  provided  smooth  surfaces 
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which  helped  minimize  flow  disturbances.  Each  wall  section  was  made 
separately  and  then  joined  to  form  the  model.  This  allowed  the  de- 
sired shape  to  be  obtained  and  also  allows  for  construction  of 
differing  geometries  in  future  experimentation. 

The  joints  between  mating  plexiglass  sections  had  to  be 
liquid  tight  and,  in  addition,  it  was  desired  that  it  be  possible 
to  disassemble  the  model.  Thus,  a permanent  chemical  bond  of  the 
plexiglass  was  not  possible.  Instead,  satisfactory  joints  were 
formed  by 

a.  Machining  all  mating  surfaces  with  an  end  mill  to  form 
flat,  square  surfaces.  This  also  resulted  in  adherence 
to  the  specified  dimensions. 

b.  Drilling  and  tapping  the  walls  and  cover  plates  to  allow 
for  joining  with  nylon  machine  screws.  Use  of  these 
screws  was  limited  somewhat  in  the  test  section  since  it 
was  desired  to  allow  for  future  changes  in  the  geometry 
of  that  section. 

c.  Both  joining  and  sealing  ail  mating  surfaces  with  silicone 
rubber  sealer. 

The  cover  plate  was  V plexiglass  plate.  Each  wall  section 
was  of  h“  plexiglass  plate  and  was  2"  in  height.  With  this  height 
and  an  inlet  width,  b,  of  V the  aspect  ratio  is  four  to  one.  Selec- 
tion of  this  aspect  ratio  and  height  was  based  on  examination  of  the 
results  presented  by  Han  [9].  According  to  Han,  for  this  aspect 
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ratio  the  flow  in  the  midplane  of  the  inlet  would  be  nearly  two- 
dimensional  . 

The  nature  of  the  cavity  flow  was  unknown  before  testing. 
Testing  has  shown  that  the  jet  flow  in  the  central  cavity  region 
remains  largely  two-dimensional.  However,  three-dimensional  effects 
were  visible  in  the  test  section  under  certain  conditions.  In  par- 
ticular, the  following  cases  were  observed. 

a.  In  the  recirculation  regions  dye  initially  on  the  mid- 
plane might  move  in  the  vertical  direction  at  points 
where  two  flow  patterns  intersected.  For  example,  such 
a point  occurred  for  low  Reynold  number  cases  where  the 
main  recirculation  flow  intersected  the  main  jet  flow. 

b.  In  the  boundary  layers  in  the  splitter  region  the  flow 
became  three  dimensional.  Although  visible  in  some 
cases,  this  boundary  layer  migration  near  the  splitter 
edge  had  little  effect  on  the  centerline  flow. 

c.  Some  combination  of  heating  from  the  lighting  for  the 
photographs  and  positioning  of  the  storage  reservoirs 
caused  a distinct  vertical  bow  in  the  main  jet  flow.  The 
exact  cause  of  this  bow  was  not  determined.  By  minimal 
use  of  the  lighting  and  by  maintaining  the  reservoirs 

at  least  several  inches  above  the  model,  the  problem 
was  avoided  altogether. 

The  base  plate  was  a 22"x51"  rectangular  plate  of  *5"  plexi- 
glass. Due  to  the  leakage  problem  it  was  necessary  to  mount  this 
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plate  on  a flat  solid  base.  As  shown  in  Fig.  3.1,  the  base  plate 
and  model  rest  on  an  angle  iron  table  frame.  Across  the  top  of  this 
frame  were  placed  a series  of  2"x4"  boards,  with  each  board  resting 
on  its  2"  edge.  These  were  planed  Mat  before  installation  and 
shimmed  on  mounting  to  produce  a flat  base.  The  plexiglass  base 
plate  was  mounted  on  this  base. 

The  model  itself  consisted  of  a screen  section,  converging 
nozzle  section,  inlet  section,  test  section,  and  two  outlet  sections. 

The  flow  entered  the  screen  section  through  two  V holes  in  the 
plexiglass  end  plate.  The  screen  section  was  12"  wide  by  63/4"  long 
and  contained  5 wire  mesh  screens,  as  shown  in  Fig.  3.3b.  Three  of 
these  screens  were  of  30  x 40  mesh  brass  screen  and  two  were  of  40x60 
mesh  brass  screen.  Each  screen  was  cold  soldered  on  a 2"  x!2"  rec- 
tangular frame  made  of  brass  brazing  rod.  Each  end  of  this  frame  was 
supported  in  a slot  machined  in  the  plexiglass  side  walls.  The  pur- 
pose of  the  screens  was  to  damp  out  disturbances  and  thus  produce  a 
very  low  turbulence  inflow. 

Disturbances  were  further  damped  by  a smooth  acceleration  I 

through  a highly  converging  nozzle,  as  shown  in  Fig.  3.2  and  Fig. 

3.3b.  This  nozzle  was  formed  by  two  sine  curves  machined  »n  plexi- 
glass blocks.  The  flow  area  changed  continuously  from  12  inihes  by 
2 inches  to  V2  inch  by  2 inches,  a 24  to  1 nozzle  ratio. 


The  converging  nozzle  was  followed  by  a V2  inch  wide  inlet 
channel.  This  section  had  a length  of  20  inches.  According  to 
Han's  calculations  [9],  this  length  would  provide  fully  developed 
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profiles  at  the  inlet  to  the  test  section.  These  inlet  profiles  were 
expected  to  be  fully  developed  for  Reynolds  numbers  up  to  approximately 
500. 

The  test  section  was  in  the  exact  shape  of  the  region  modeled 
in  the  computer  simulation  and  is  visible  in  Fig.  3.3a.  The  cavity 
region  was  a 5 inch  square.  The  receiver  outlets  were  rounded  to  pro- 
vide a smooth  transition  into  a relatively  large  outlet  section  (423/St 
inches  by  5 inches).  This  configuration  minimized  the  downstream 
effects  on  the  receiver  outlet  flow.  The  flow  left  each  outlet  sec- 
tion through  a 3A  inch  diameter  hole  positioned  in  the  center  of  the 
wall  opposite  each  receiver  outlet. 

Provisions  were  made  on  the  model  for  the  mounting  of  dye  in- 
jection fittings,  which  will  be  described  later,  and  air  bleed  fittings. 
Except  as  noted  later,  these  fittings  were  mounted  in  holes  drilled  in 
the  cover  plate.  The  bleed  holes  were  positioned  at  corner  points  in 
the  test  and  outlet  sections  and  were  used  to  remove  air  trapped  under 
the  cover  plate  when  filling  the  model  with  fluid. 

C.  Fluid  Supply  System 

Flow  to  and  from  the  model  was  controlled  through  movement  of 
a system  of  reservoirs.  The  actual  flow  through  the  model  depends 
on  the  height,  or  pressure,  differential  between  the  feed  reservoir 
and  the  two  outlet  reservoirs.  This  pressure  differential  was  changed 
by  vertical  movement  of  one  or  more  of  the  reservoirs.  The  overall 
arrangement  is  shown  in  Fig.  3.3c. 
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The  feed  reservoir,  the  larger  reservoir  on  the  left  in  Fig. 

3.3c,  had  a second  purpose.  It  served  to  shield  the  flow  to  the 

model  from  disturbances  created  by  the  pump  which  forced  the  fluid 

to  the  level  of  the  feed  reservoir.  The  feed  reservoir  acted  as  a 
stilling  chamber. 

All  three  reservoirs  were  rectangular  boxes  made  of  V2  inch 
and  % inch  plexiglass.  The  plexiglass  sides  and  base  were  milled 
and  fastened  together  with  nylon  machine  screws.  All  joints  were 
then  sealed  through  use  of  ethylene  dichloride  to  form  a chemical 
bond  between  mating  plexiglass  ports.  All  three  reservoirs  were 
divided  into  two  compartments  by  a vertical  wall,  or  wier.  One  com- 
partment in  each  reservoir  drained  back  to  the  storage  reservoir  and 
on  typical  runs  was  nearly  empty. 

In  the  feed  reservoir  the  main  compartment  had  two  ports. 

The  flow  from  the  supply  pump  entered  through  one  port  and  left 
through  the  other  port  to  flow  to  the  model.  Excess  fluid  from  the 
supply  pump  flowed  over  the  wier  and  back  to  the  storage  reservoir. 
The  head,  or  level,  in  the  reservoir  was  thus  fixed  at  the  level  of 
the  top  of  the  wier. 

The  reservoirs  were  supported  by  the  guide  bracket,  cable, 
and  pipe  arrangement  shown  in  Fig.  3.3c.  The  pipe  frame  was  of  1 V4 
inch  diameter  standard  steel  pipe.  One  end  of  each  pipe  was  flanged 
to  the  floor  and  the  other  end  is  flanged  to  a section  of  4 inch 
steel  channel.  The  frame  was  further  supported  by  sections  of  3A 
inch  pipe  which  connect  the  vertical  1 Vi,  inch  pipe  to  the  room  wall. 


The  guide  brackets  slid  vertically  on  the  1 ’/■,  inch  pipe  and 
served  to  position  the  reservoirs  in  the  horizontal  plane.  The  weight 
of  the  reservoirs  and  guide  brackets  was  supported  by  3/K  inch  plastic 
coated  tiller  cable,  which  is  visible  in  Fig.  3.1.  Each  of  these 
cables  ran  from  the  guide  bracket  and  reservoir  through  the  upper 
channel,  over  two  rollers  mounted  in  the  channel,  and  down  the  other 
side  of  the  1 'A  inch  pipe.  At  the  floor  the  cable  makes  a 180°  turn 
through  a cable  pulley  and  was  clamped  to  itself  with  a cable  clamp. 
Loosening  of  this  clamp  allowed  coarse  adjustment  of  reservoir  height. 
Fine  adjustment  was  provided  by  a turnbuckle  which  was  mounted  between 
the  tiller  cable  and  the  reservoir  guide  bracket.  One  full  turn  of 
this  turnbuckle  moved  the  reservoir  vertically  approximately  0.1 
inches. 

Flow  between  the  reservoirs,  model,  pump,  and  storage  reser- 
voirs is  through  1 inch  Gates  heater  hose  where  the  line  must  be 
flexible  and  through  1 inch  PVC  pipe  where  a solid  line  was  accept- 
able. The  fluid  was  pumped  from  the  storage  reservoir  to  the  feed 
reservoir  by  a Model  DH11  Eastern  Industries  centrifugal  pump. 

D . Instrumentation 

The  measurement  of  the  flow  rate  was  complicated  by  the  small 
flow  rates  involved.  It  was  decided  to  construct  a flow  section 
across  which  the  pressure  drop  could  be  easily  measured  and  corre- 
lated with  the  flow  rate.  The  resulting  design  is  shown  in  Fig.  3.3d. 
The  design  consisted  or  two  brass  nozzles  connected  by  a length  of  Vi* 
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inch  I.D.  tygon  tubing.  The  tubing  length  was  set  to  yield  a suitable 
pressure  drop  for  the  flow  range  expected.  The  pressure  drop  was 
measured  on  a manometer  board  having  a range  of  1 to  12  inches  of 
water. 

Calibration  for  a given  fluid  and  length  of  tubing  was  accom- 
plished easily.  The  meter  was  connected  in  a line  between  the  feed 
reservoir  and  one  of  the  outlet  reservoirs.  These  reservoirs  were 
then  adjusted  to  produce  different  flow  rates.  Each  flow  rate  pro- 
duced a different  pressure  drop,  which  was  recorded  along  with  the 
flow  rate.  To  determine  the  flow  rate,  the  flow  from  the  outlet  res- 
ervoir to  the  storage  reservoir  was  diverted  for  a measured  period  of 
time  and  collected  in  a pail.  The  collected  fluid  was  weighed  on  a 
balance  beam  scale,  which  measured  in  increments  of  0.01  lb.  The 
time  period  was  measured  with  a stopwatch,  which  measured  in  incre- 
ments of  0.01  minutes.  The  weight  flow  rate  was  then  calculated  by 
dividing  the  fluid  weight  by  the  measured  time.  Typical  runs  were 
for  five  to  ten  minutes  and  collected  between  2 and  6 pounds  of 
fluid.  Calibration  curves  for  the  flow  meters  used  are  shown  on 
Fig.  3.4.  These  curves  are  for  water  at  approximately  73°F.  During 
the  experiments  the  temperature  varied  slightly  from  this  temperature 
(±3°F),  but  the  effect  of  this  variation  on  the  flow  rates  measured 
would  be  small  and  was  ignored. 

Two  flow  meters  were  used  to  obtain  the  inlet  and  one  of  the 
outlet  flows.  The  other  outlet  flow  was  taken  as  the  difference 


between  these  two  measured  flows. 


Figure  3.4.  Flow  Meter  Calibration  Curves 
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Once  the  flow  rate  was  determined  it  had  to  be  related  to  the 
inlet  Reynolds  number  used  in  the  computer  solutions.  The  average 
velocity  in  the  inlet  channel,  is  given  by  the  volumetric 

flow  rate,  Q,  divided  by  the  channel  area,  A . From  Han's  calcula- 
tions [9],  the  centerline  velocity,  Uq,  for  a channel  with  an 
aspect  ratio  of  four  is 


— = 1.7728  . (3.1a) 

3-D 


At  the  horizontal  midplane  the  velocity  profile  is  parabolic  with  a 
maximum  velocity  of  uQ.  For  the  computer  solutions  the  velocity  on 
which  the  Reynolds  number  was  based  was  the  average  velocity  of  a 
two  dimensional  parabolic  profile  for  which 


u0  ' 1,5  U2-D 


(3.1b) 


Therefore,  the  velocity  for  use  in  the  Reynolds  number  is 
U2-D  T75  1-182  U3-D  • 


(3.1c) 


Then , 


U2-0  * ,J82^ 


and 


Re 


1.182  Q 
h v 


(3. Id) 


(3.  le) 


where  h is  the  height  of  the  inlet  channel. 
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As  mentioned  earlier  the  flow  visualization  was  accomplished 
by  injecting  dye  into  the  flow.  Initially  there  was  some  concern 
that  the  needles  used  for  injecting  the  dye  would  disturb  the  flow 
patterns.  For  this  reason  the  hydrogen  bubble  flow  visualization 
technique  [26]  was  investigated.  This  technique  involves  generating 
hydrogen  bubbles  by  the  hydrolysis  of  water  on  a small  wire  placed 
in  the  flow.  The  size  of  the  wire  is  such  that  the  disturbance  to 
the  flow  would  be  minimal.  The  hydrogen  bubbles  are  swept  off  the 
wire  and  heir  movement  indicates  the  flow  patterns.  Unfortunately 
it  was  found  that  for  the  low  flow  rates  considered  in  this  study 
the  bubbles  would  rise  and  collect  beneath  the  cover  plate.  It  was 
expected  that  the  bubbles  would  go  back  into  solution,  but  this  was 
not  the  case,  and  the  hydrogen  bubble  technique  was  abandoned.  Dye 
injection  was  then  examined  and  it  was  found  that  the  injector  did 
not  noticeably  disturb  the  flow  field. 

The  fluid  used  for  the  flow  visualization  was  water.  The 
temperature  was  measured  using  a standard  mercury  thermometer.  The 
fluid  properties  were  then  taken  as  those  of  pure  water  at  this 
temperature.  This  introduces  some  error  into  the  Reynolds  number 
determination  because  the  water  is  not  pure.  However,  the  desired 
results  are  qualitative  and  it  was  felt  that  this  error  is  not  sig- 
ni ficant. 

The  dye  chosen  was  a mixture  of  red  and  green  food  coloring 
and  water,  proportioned  to  produce  a dark  navy  blue  or  black  color. 
The  dye  was  mixed  at  the  time  of  each  experiment  with  water  taken 





from  the  storage  reservoir.  This  insured  that  the  dye  and  the  water 
in  the  flow  were  at  the  same  temperature  and  thus  would  have  approxi- 
mately the  same  density. 

The  dye  was  introduced  into  the  flow  through  two  separate  in- 
jectors. In  one  injector  the  dye  was  injected  through  a 0.032  inch 
O.D.  by  0.020  inch  I.D.  tube  which  was  held  in  place  by  soldering  it 
into  a hose  fitting.  The  tube,  which  was  mounted  vertically,  had  one 
end  bent  in  order  that  the  dye  would  be  injected  horizontally  and  at 
the  model  midplane.  This  injector  was  mounted  such  that  dye  was 
introduced  into  the  inlet  channel  approximately  3 inches  from  the 
test  section. 

The  second  injector  was  similar  in  construction  to  the  first 
except  that  the  tube  was  not  bent.  This  injector  was  mounted  in  a 
side  wall  of  the  test  section  such  that  the  tube  lay  completely  in 
the  horizontal  midplane.  The  tube  was  held  in  place  in  such  a way 
that  it  could  be  easily  moved  in  a direction  normal  to  the  wall. 

This  allowed  visualization  of  recirculation  streamlines  over  the 
entire  recirculation  region. 

For  both  injectors  a Vi6  inch  I.D.  tygon  tube  was  attached 
to  the  injector  and  extended  to  a fitting  which  was  soldered  to  a 
hypodermic  needle.  The  hypodermic  needle  was  placed  on  a syringe 
which  contained  the  dye.  The  dye  could  then  be  injected  into  the 
model  in  two  ways.  It  could  be  injected  by  simply  forcing  the 
syringe  plunger  by  hand.  A1 ternati vely , the  syringe  could  be  ele- 
vated, say  by  supporting  it  with  a ringstand,  and  for  the  proper 
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elevation  the  dye  would  be  forced  into  the  flow  through  a combina- 
tion of  the  weight  of  the  dye  and  of  the  syringe  plunger. 

E . Photographic  Equipment 

The  streamlines  shown  by  the  dye  patterns  were  recorded  on 
film.  The  main  considerations  were  the  choice  of  a camera  and  the 
method  of  1 ighting. 

The  camera  chosen  was  a Canon  Canonet  QL  17.  This  model  is 
a single  lens,  35  mm  camera,  with  this  particular  one  having  a 1:1.7 
lens.  A +3  close-up  lens  was  used  in  conjunction  with  the  camera. 

The  camera  was  mounted  above  the  model  on  a horizontal 
support  fixed  to  one  of  the  pipes  in  the  reservoir  support  frame. 

The  support  was  an  aluminum  tube  which  formed  a cantilever  beam 
extending  over  the  model.  The  height  of  the  beam  above  the  model 
was  controlled  by  the  area  desired  in  the  picture,  the  close-up  lens 
used,  and  the  footage  setting  on  the  camera.  For  the  +3  lens,  a 
footage  setting  of  15  feet,  and  a desired  coverage  of  a 10"  by  7" 
area,  the  height  from  model  midplane  to  lens  was  12 'A  inches. 

The  lighting  was  furnished  by  a single  500  watt  photoflood 
bulb  mounted  in  a standard  circular  reflector.  This  was  centered 
directly  below  the  test  section  and  positioned  approximately  28 
inches  below  the  test  section  midplane.  This  strong  beam  of  light 
was  diffused  by  passage  through  a diffuser  made  of  vellum  paper. 

The  diffuser  was  positioned  4 inches  below  the  model. 
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F . Comments  on  Experimental  Procedure 

Use  of  the  above  equipment,  excluding  the  photographic  equip- 
ment, was  rather  straightforward.  The  photographic  equipment  also  was 
not  complex  as  long  as  standard  directions  for  the  respective  pieces 
of  equipment  were  available  and  followed.  The  general  procedure 
for  use  of  each  piece  of  equipment  should  be  apparent  from  the 
description  of  the  equipment. 

There  were,  however,  several  experimental  techniques  which  are 
worthy  of  mention.  These  were  as  follows: 

a.  Care  had  to  be  exercised  to  insure  that  all  the  air  had 
been  removed  from  the  model  and  the  supply  lines.  This 
problem  was  alleviated  by  inclining  the  model  as  it 
filled  with  water.  This  allowed  the  air  to  rise  to  the 
bleed  holes.  Trapped  air  could  be  eliminated  from  the 
supply  lines  by  making  sure  that  the  high  points  in 
these  lines  were  points  at  which  air  could  be  bled. 

b.  The  feed  and  outlet  reservoirs  interacted  in  their  effect 
on  the  flow  rate.  For  example,  blocking  one  receiver 
would  not  result  in  all  the  previous  flow  going  through 
the  other  receiver.  Instead,  the  receiver  flow  rate 
would  increase  and  the  inlet  flow  rate  would  decrease. 
Thus,  some  manipulating  was  necessary  to  obtain  the  de- 
sired flow  rates  when  changing  from  one  combination  of 
flow  rates  to  another. 
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c.  The  time  constants  of  this  system  were  rather  long. 

Thus,  after  changes  in  flow  rates,  a sufficient  period 
of  time  had  to  be  allowed  for  the  system  to  come  to 
equilibrium.  For  low  Reynolds  numbers  this  time  period 
was  commonly  15  minutes. 

d.  The  injection  of  dye  was  somewhat  of  an  art.  Once  a 
flow  pattern  was  established,  several  attempts  at  obtain- 
ing a suitable  dye  pattern  might  be  required.  Also  the 
dye  patterns  were  not  extremely  stable  and  could  easily 
be  effected  by  vibration  of  the  model  or  even  over- 
heating from  the  lighting  for  the  camera,  as  mentioned 
earl ier. 

G.  Sample  Result 

A typical  result  for  the  dye  injection  studies  is  shown  in 
Fig.  3.5.  Here  the  flow  had  a Reynolds  number  of  100,  which  corre- 
sponds to  an  average  inflow  velocity  of  0.3  inch  per  second.  The 
dye  was  injected  to  illustrate  the  center  streamline  for  a symmetric 


case. 


CHAPTER  IV 


DISCUSSION  OF  RESULTS 

Tn*1  I..  ■.  Leal  node!  developed  in  Chapter  II  was  used  to  pre- 
dict the  flow  fas  and  pressure  distributions  for  a variety  of  flow 

situation-.  -Lions  were  made  in  Reynolds  number,  receiver  flows, 

and  deflect1-'  of  the  inlet  jet.  The  resulting  flow  and  pressure 
distribution  1 be  presented  and  discussed  in  this  chapter.  In 
addition,  s cr-nf  ine  photographs  were  obtained  using  the  experimental 
apparatus  d-  . ed  n Chapter  III  for  three  of  the  cases  including 
asymmet  !rn  ...q.  Ihese  photographs  will  be  compared  to  the  corre- 
sponding nun:*  a reamline  predications. 

one  of  the  projected  uses  of  the  numerical  model  is  as  a de- 
sign too  if  obtaim .ig  characteristic  curves,  as  mentioned  in  Chapter 
I.  Gene,  it  ■ o.  such  a curve  is  described  and  an  example  presented. 
This  is  f c I owed  discussion  of  the  experimental  observation  of 

The  chapter  concludes  with  several  comments 
on  nut*  cal  . aracteristics  of  the  computer  model  itself. 

A.  Flow  Fi  ; pits  for  the  Undeflected  Jet 

The  predicted  flow  fields  are  conveniently  presented  as  con- 
tour plo*  team  function  values.  Such  a plot  is  shown  in 

Fig.  4.1.  Each  solid  line  in  Fig.  4.1  represents  a line  of  constant 
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ijj , i.e.  a streamline.  Each  streamline  divides  the  flow  into  separate 
regions  since  flow  is  tangent  to  a streamline.  Thus,  the  0 and  -1.0 
streamlines  bound  the  inflow  from  the  inlet  nozzle  and  represent 
lines  of  separation  between  the  flow  entering  the  cavity  and  the  flow 
which  remains  in  the  recirculation  regions.  In  the  main  flow  the 
streamlines  vary  by  increments  of  0.125.  In  the  recirculating  flow 
the  streamlines  are  labeled,  with  the  corresponding  values  for  ip 
presented  in  Table  4.1. 

As  denoted  by  the  arrow,  the  flow  enters  the  amplifier 
through  the  single  port  in  the  lower  boundary.  For  the  flow  shown, 
which  has  a Reynolds  number  of  1,  the  jet  spreads  as  it  crosses  the 
cavity  and,  due  to  the  strong  diffusion  effect  at  this  Reynolds  num- 
ber, occupies  most  of  the  cavity  region.  The  flow  is  almost  stagnant 
in  the  upper  cavity  corners  and  in  the  lower  corners  the  jet  drives 
two  small  vorticies,  or  recirculation  regions.  The  velocities  are 
relatively  large  in  the  central  cavity  region,  where  the  x-gradient 
of  ip  is  relatively  large.  The  velocities  are  much  lower  on  the  edges 
of  the  main  flow  and  in  the  recirculation  regions.  In  both  cases 
this  is  indicated  on  Fig.  4.1  by  the  small  spatial  gradients  in  \p  in 
these  regions.  In  Fig.  4.1  one  half  of  the  flow  exits  through  each 
receiver.  Since  the  flow  is  initially  undeflected,  the  stream  func- 
tion distribution  is  symmetric. 

A series  of  such  symmetric  flows  is  shown  in  Fig.  4.2.  In 
this  series  the  Reynolds  numbers  are  respectively  20,  100,  275,  and 
1000.  Increasing  the  Reynolds  number,  Ub/v,  corresponds  to 
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• t velocity,  U,  or  decreasing  the  kinematic 
r Utions  imply  that  the  importance  of  dif- 
• is  reduced  and  that  the  flow  becomes  domi- 
ter-fiis  in  Eq.(2.6).  As  pointed  out  for  Re  = 1 , 
reads  to  the  width  of  the  cavity  region, 
ig  4.2,  the  spread  of  the  main  flow  becomes 

,'i'i  How  has  sufficient  momentum  to  drive  one 
lion  in  each  half  of  the  cavity.  For  Re  = 100 
’eci rculation  vortex  become  larger  and  the  re- 
rv.  to  drive  secondary  vorticies  in  the  cavity 
a, id  4. 2d  the  higher  momentum  in  the  main 
d:  -e  progressively  larger  secondary  vortic- 
; what  'light  be  physically  expected.  In  other 
item  of  the  entire  region  is  increased  and  the 
•odicate  the  redistribution  of  this  momentum 
-c  lating  vorticies. 

eo  in  Section  F,  the  flow  pattern  for 
-i  entally  obtainable  due  to  edgetone  oscilla- 
t.  show  that  the  numerical  method  was  stable 
■•e  / t i:  value  for  Reynolds  number.  The  upper  bound 
i afore  the  onset  of  oscillation  with  the  partic- 
d s used  in  this  study  is  roughly  275. 
l hod  of  presenting  the  flow  field  results  is 
' on  for  different  sections  in  the 


showr  Fig.  4.3a.  This  is  a v-velocity 
- ; is  normal  to  the  amplifier  centerline  at 

; >3;  which  corresponds  to  a y-di stance 
t • ■ • . o'cttd  for  three  Reynolds  numbers 
symmetry,  flows.  Again  the  effect  of  in- 
’s 1 early  visible.  For  the  higher  Reynolds 
, rile  h > narrowed  and  a higher  centerline 

The  momentum  of  the  inlet  flow  has  not 
i ' ir  o the  sides  of  the  cavity  region. 

. . : ol  diffusion  is  visible  in  Fig.  4.3b.  The 
.•  ■ )cities  for  1 < i < 41  and  j = 39,  which 

o inyes  on  the  splitter  and  is  located  at  a 
For  the  low  Reynolds  number  profile  the  presence 
ii.i  w,  stream,  a diffusion  phenomenon,  and 
. 1 e shown  with  the  centerline  velocity  being 
not  observed  fur  the  higher  Reynolds  number 
■ . i:.  reduction  is  observed,  the  domination  of 

<e  ion  prevents  any  distortion  of  the  profile, 
ofile  plotted  on  Fig.  4.3a  is  that  for  a laminar 
. ofile  might  be  used  in  conjunction  with  the 
1].  The  profile  is  based  on  Scnl ichting ’ s 
jer.  [25]  The  constants  have  been  evalu- 
. piv; i 7 e is  non-dimensional.  It  also  has  the  same 

1 at  the  amplifier  inlet  as  the  numerical 
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solutions.  The  laminar  jet  profile  is  for  a Reynolds  number  of  100 
and  thus  can  be  compared  to  tne  Re=100  velocity  profiles  from  the 
numerical  results. 

It  was  not  expected  that  the  profiles  would  match,  since  it 
is  assumed  in  the  derivation  of  the  laminar  jet  profile  that  the  jet 
exhausts  into  an  infinite  medium  and  that  its  momentum  is  constant. 
However,  the  general  shapes  are  similar,  with  the  free  jet  having  a 
lower  velocity  at  all  points  in  the  central  flow  region.  Thus,  for 
the  central  region  the  momentum  of  the  amplifier  flow  did  not  remain 
constant;  it  increased.  This  indicates  that  the  positive  v-velocity 
component  of  the  recirculation  flow  tends  to  increase  the  momentum 
of  the  main  flow,  even  though  the  main  flow  is  driving  the  recircula- 
tion flow.  Thus,  the  confinement  of  the  jet  by  the  amplifier  walls 
concentrates  the  positive  y-momentum  flux  in  the  central  flow  region. 

A further  comparison  of  velocity  profiles  is  shown  in  Fig. 
4.4.  This  plot  illustrates  the  centerline  velocity  decay  for  the 
three  Reynolds  numbers  and  for  the  free  laminar  jet  with  Re =100. 

As  the  Reynolds  number  is  increased  the  viscous  effects  on  the  main 
flow  are  reduced,  which  results  in  the  centerline  velocity  remaining 
high.  Again,  the  laminar  free  jet  velocity  decays  more  rapidly  than 
the  corresponding  amplifier  flow  velocity  profile  due  to  the  effect 
the  recirculation  regions  have  on  the  main  flow. 

As  discussed  in  Chapter  III,  an  effort  was  made  to  obtain 
experimental  verification  of  the  predicted  flow  fields.  Flow  field 
photographs  were  obtained  for  two  of  the  symmetric  cases,  Re = 100 
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and  Re =20,  and  for  one  unsymmetric  case.  The  latter  case  will  be 
presented  later  in  this  chapter. 

A photograph  of  the  flow  *or  Re =100  is  presented  in  Fig.  4.5 
along  with  the  corresponding  numerically  predicted  streamline  pattern. 
This  visualization  was  obtained  by  injecting  a large  amount  of  dye 
into  the  cavity  region,  disrupting  the  flow  pattern,  and  allowing 
the  flow  pattern  to  redevelop.  Except  for  the  large  concentration 
of  dye  located  on  the  base  plate  near  the  splitter  this  technique 
worked  well.  All  the  dye  is  certainly  not  on  the  midplane,  but  the 
flow  pattern  is  evident.  The  dye  concentration  mentioned  above  had 
settled  into  the  lower  boundary  layer  and  had  not  had  time  to  wash 
out  when  the  picture  was  taken.  The  center  streamline  is  on  the  mid- 
plane and  slightly  off  center  horizontally.  It  is  a continuous 
stream  being  injected  through  a needle,  the  injection  technique 
described  in  Chapter  III. 

A comparison  between  (a)  and  (b)  of  Fig.  4.5  is  felt  to  be 
excellent  qualitative  verification  of  the  computer  flow  solution. 

In  Fig.  4.5b  the  boundaries  of  the  main  flow  are  clearly  indicated, 
as  is  the  flow  pattern  in  the  main  recirculation  region. 

Details  of  the  flow  shown  in  Fig.  4.5b  were  obtained  on  the 
photographs  shown  in  Fig.  4.6.  Again,  these  compare  well  with  the 
stream  function  plot  in  Fig.  4.5a.  For  Fig.  4.6  dye  was  injected 
into  the  cavity  midplane  through  the  needle  which  was  extended  from 
the  left  wall.  This  needle  was  repositioned  between  each  picture  to 
produce  the  different  streamlines  shown. 
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a)  Central  Vortex  Region 
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Figure  4.6a  shows  a central  vortex  streamline  which  completes 
a closed  vortex  loop.  Figure  4.6b  shows  another  recirculation  stream- 
line which  is  farther  from  the  center  of  the  vortex.  In.  Fig.  4.6c 
the  streamline  shown  is  in  a region  of  very  slow  flow  and  most  of  the 
main  recirculation  flow  moves  between  the  position  of  this  streamline 
and  the  center  of  the  vortex.  The  streamline  shown  in  Fig.  4.6d 
marks  the  point  where  the  recirculation  pattern  separates  from  the 
wall.  The  flow  illustrated  here  was  extremely  slow  and  the  stream- 
line would  not  extend  farther  into  the  flow  field  without  falling 
into  the  lower  boundary  layer  and/or  diffusing.  This  is  consistent 
with  the  low  velocities  indicated  in  Fig.  4.5a.  Further  details  on 
this  flow  field  will  be  presented  later  in  conjunction  with  compari- 
sons between  symmetric  and  unsymmetric  flows  for  Re =100. 

A comparison  similar  to  that  of  Fig.  4.5  is  made  in  Fig.  4.7 
for  the  symmetric  case  with  Re = 20.  For  the  experimental  model  used, 
having  b = 0. 5 inches,  the  flow  rate  at  this  Reynolds  number  was  ex- 
tremely low,  0.12  lb  /min.  For  reference,  at  this  flow  rate  it 
takes  70  minutes  to  fill  a one  gallon  container.  The  low  velocities 
of  the  flow  made  the  visualizations  difficult  to  obtain. 

In  Fig.  4.7b  the  flow  pattern  in  the  right  cavity  was  ob- 
tained by  flooding  this  portion  of  the  cavity  with  dye  and  then 
allowing  the  flow  pattern  to  redevelop.  In  the  left  cavt. 
an  individual  streamline  is  shown.  Due  to  the  low  ve  • • 
slight  density  differences  involved,  this  str*  • • I un- 

toward the  base  plate  after  traveling  th« 
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Although  the  streamlines  in  Fig.  4.7b  are  not  as  well  defined 
as  in  Figs.  4.5b  and  4.6,  several  observations  about  the  flow  field 
can  be  made.  First,  in  comparing  Figs.  4.5b  and  4.7b,  it  is  observed 
that  the  center  of  the  vortex  region  has  moved  towards  the  inlet. 

This  movement  is  in  agreement  with  the  contour  plots  shown  in  Fig. 
4.5a  and  Fig.  4.7a.  Second,  the  recirculation  pattern  in  the  right 
half  of  the  cavity  region  in  Fig.  4.7b  is  in  qualitative  agreement 
with  the  contour  plot.  In  particular,  the  experimental  vortex  center 
appears  to  be  located  in  the  position  predicted.  Third,  the  boundary 
between  the  main  flow  and  the  recirculation  region  is  in  qualitative 
agreement.  In  comparison  to  the  Re =100  case,  the  jet  boundary  in 
the  right  half  cavity  has  moved  to  the  right.  This  implies  a larger 
spread  of  the  main  flow:  a result  of  the  relative  increase  in  diffu- 
sion at  Re  = 20. 

As  discussed  in  Chapter  II,  provisions  were  made  in  the  nu- 
merical simulation  for  varying  the  flow  rates  in  the  receivers  while 
keeping  the  inflow,  or  Reynolds  number,  constant.  Such  a series  with 
Re =100  is  shown  in  Fig.  4.8.  Here  the  flow  rate  in  the  left  re- 
ceiver is  0.5  Q$,  0.625  Qs , 0.750  Qg,  and  1.0  Q$.  The  dominant  fea- 
ture of  these  plots  is  a progressively  larger  shift  of  the  main  flow 
boundaries  to  the  left.  This  would  be  expected  because  progressively 
more  flow  is  moving  through  the  left  receiver.  The  movement  of  the 
main  flow  is  accompanied  by  a shift  of  the  center  of  the  main  recir- 
culation regions  to  the  left.  Also,  the  total  mass  circulating  in 
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the  left  vortex  region  is  reduced.  This  is  indicated  by  the  reduc- 
tion in  the  maximum  value  of  the  stream  function  found  in  the  left 
vortex  as  the  series  progresses  from  Fig.  4.8a  to  4.8d. 

Photographs  were  obtained  for  the  blocked  case,  i.e.  QL  = 

1.0  Q$,  in  this  series.  The  first  of  these  photographs  is  shown  in 
Fig.  4.9,  along  with  the  corresponding  contour  plot.  This  visualiza- 
tion was  made  by  flooding  the  region  with  dye  and  waiting  until  the 
flow  patterns  had  redeveloped.  In  addition,  more  detail  was  obtained 
by  injecting  from  the  left  wall  while  the  flow  pattern  was  redevelop- 
ing in  the  right  cavity  region.  This  resulted  in  a poorly  defined 
pattern  in  the  central  vortex  region  and  a well  defined  streamline 
along  the  lower  edge  of  the  main  recirculation  region. 

Several  observations  can  be  made  from  this  comparison.  First, 
the  boundary  between  the  right  recirculation  region  and  the  main  flow 
is  well  defined  and  its  location  agrees  well  with  the  numerical  re- 
sults shown  in  the  contour  plot.  Second,  the  centers  of  the  recir- 
culation regions  can  be  inferred  and,  within  the  limits  of  such  an 
inference,  agree  with  the  center  locations  predicted.  Third,  the 
well  defined  individual  streamline  in  the  left  cavity  is  almost  iden- 
tical to  one  of  the  streamlines  indicated  on  the  contour  plot  for  the 
same  region.  Fourth,  near  the  main  flow  boundary  in  the  lower  right 
recirculation  region  the  shape  of  the  streamlines  can  be  visualized. 
Again  this  shape  compares  favorably  with  the  predicted  flow  pattern. 
Finally,  a flow  into  the  right  receiver  is  faintly  indicated  on  the 
photograph.  However,  it  is  not  certain  whether  this  flow  was  on  the 


86 

midplane  or  in  the  boundary  layer  on  the  base  plate.  If  on  the  mid- 
plane it  would  be  in  agreement  with  the  predicted  flow. 

Two  photographs  of  flow  details  for  the  blocked  case  are 
shown  in  Fig.  4.10b  and  Fig.  4 . 1 Od . Both  are  streamlines  in  the 
main  flow.  One,  Fig.  4.10b,  shows  the  center  inlet  streamline  and 
the  other.  Fig.  4 . 1 Od , shows  a streamline  very  near  the  inlet  wall. 

Both  agree  qualitatively,  except  as  noted  shortly,  with  the  corre- 
sponding numerical  streamlines  shown  in  Fig.  4.9a. 

Although  it  would  be  desirable,  it  is  difficult  to  get  quan- 
titative comparisons  between  prediction  and  experiment  due  to  both 
paralax  problems  and  measurement  problems.  Therefore,  this  was  not 
attempted  to  any  extent  in  this  study.  However,  an  approximate 
quantitative  comparison  can  be  obtained  from  Fig.  4 . 1 Od  because  the 
camera  is  placed  almost  directly  over  the  splitter  tip.  Measurement 
from  the  photograph  indicates  that  the  streamline  enters  the  receiver 
at  a point  located  roughly  61%  of  the  receiver  width  from  the  split- 
ter tip.  Measurement  on  the  contour  plot  locates  the  same  point  at 
a distance  of  approximately  66%  of  the  receiver  width. 

Close  examination  of  the  streamline  in  Fig.  4 . 1 Od  shows  some 
deviation  from  the  computer  plot.  In  particular,  the  computer  model 
predicts  that  the  wall  streamline,  the  zero  streamline,  is  slightly 
to  the  left  of  the  streamline  in  the  visualization.  In  fact,  in 
the  region  near  the  receivers  the  visualization  streamline  agrees 
much  better  with  the  predicted  = -0.125  streamline.  For  example. 


b)  Center  Streamline 
Blocked  Flow 


a)  Center  Streamline 
Symmetric  Flow 


d)  Wall  Streamline 
Blocked  Flow 


c)  Wall  Streamline, 

Symmetric  Flow 

Figure  4.10.  Comparison  Between  Symmetric  and 
Blocked  Flows  with  RE  = 100 
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This  deviation  can  be  partially  explained  because  the  dye  forming  the 
streamline  had  some  width  as  it  left  the  inlet;  that  is,  it  was  not 
all  right  against  the  wall.  Therefore,  it  is  not  exactly  the  zero 
streamline  and  would  not  be  expected  to  agree  exactly  with  the  com- 
puter prediction  for  the  zero  streamline.  However,  there  does 
appear  to  be  some  other  discrepancy  and  this  could  possibly  be  an 
indication  of  the  numerical  problem  to  be  discussed  in  Section  D of 
this  chapter. 

The  remaining  pictures  in  Fig.  4.10  are  the  details  of  the 
symmetric  case  with  Re  = 100  which  were  mentioned  earlier.  The  two 
details  shown  are  the  center  streamline.  Fig.  4.10a,  and  the  wall 

'i 

streamline.  Fig.  4.10c.  Two  comparisons  can  be  made.  First,  these 
streamlines  can  be  compared  with  the  corresponding  experimental 
streamlines  for  the  blocked  case;  that  is,  compare  Fig.  4.10a  to 
Fig.  4.10b  and  Fig.  4.10c  to  Fig.  4.10d.  This  clearly  indicates  the 
shift  of  the  main  flow  to  the  left  for  the  blocked  case.  Comparing 
Fig.  4.10c  and  Fig.  4 . 1 Od  also  indicates  that  the  flow  becomes  more 
complex  for  the  blocked  case.  In  the  symmetric  case  the  wall  stream- 
line enters  the  receiver  smoothly,  with  only  a slight  bend,  while  in 
the  blocked  case  a pronounced  bend  occurs  just  before  the  streamline 
enters  the  receiver. 

The  second  comparison  which  can  be  made  is  to  compare  Fig. 

4.10a  and  Fig.  4.10c  to  the  corresponding  contour  plot.  Fig.  4.5a. 

The  center  streamline  and  the  wall  streamline  agree  fairly  well. 

After  entering  the  receiver  the  wall  streamline  does  not  correspond 


exactly  with  the  receiver  wall,  but  this  might  be  expected  because, 
again,  the  experimental  streamline  was  not  exactly  the  zero  stream- 
line due  to  its  finite  width  at  the  inlet.  The  experimental  stream- 
line also  does  not  show  the  same  curvature  the  numerical  zero 
streamline  shows  as  it  terminates  on  the  receiver  wall.  Since  the 
experimental  streamline  is  not  exactly  the  zero  streamline  this 
behavior  appears  to  correspond  acceptably  with  that  predicted  numer- 
ical ly. 

B.  Flow  Field  Results  for  the  Deflected  Jet 

As  discussed  in  Chapter  II,  the  initial  deflection  of  the 
inlet  flow  was  modeled  as  an  initial  angular  displacement  $ of  the 
inlet  velocity  vector.  Sample  results  for  deflected  flows  are  shown 
in  Fig.  4.11.  The  flow  modeled  in  Fig.  4.11a  has  a Reynolds  number 
of  100,  equal  receiver  flow  rates,  and  a 5°  deflection.  Some  portions 
of  the  flow  are  noticeably  unsymmetric,  particularly  the  lower  sec- 
ondary vorticies  and  the  stream  function  magnitudes  reached  in  the 
main  vorticies.  However,  the  effect  on  the  main  flow  is  small  and 
the  flow  is  almost  symmetric  in  the  receiver  region. 

In  Fig.  4.11b  the  flow  conditions  are  the  same  as  in  Fig. 

4.11a  except  the  initial  deflection  is  increased  to  15°.  Compared 
to  Fig.  4.11a  the  effect  of  the  additional  deflection  is  pronounced. 
The  main  flow  is  noticeably  deflected  throughout  the  cavity.  Ini- 
tially there  is  a strong  deflection  to  the  left  and  then  the  flow 
bends  back  to  discharge  through  the  receivers  equally.  The 
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secondary  vorticies  in  the  left  cavity  region  have  disappeared  and 
the  unsyimetric  character  of  the  main  recirculation  region  has  in- 
creased. 

Figure  4.11c  is  presented  as  an  example  of  a deflected  flow 
with  unequal  flow  rates  in  the  receivers.  Here  the  deflection  is 
15°,  Re =100,  and  Q-j  = 0.375  Q . Compared  to  Fig.  4.11b  the  effect 
was  to  increase  the  curvature  of  the  streamlines  in  the  main  flow, 
since  they  must  bend  back  farther  to  allow  0.625  of  the  flow  to  exit 
through  the  right  receiver. 

C.  Pressure  Field  Results 

Presentation  of  the  pressure  results  divides  conveniently 
into  three  parts.  First,  a series  of  pressure  solutions  correspond- 
ing to  the  flow  cases  discussed  above  will  be  presented.  This  is 
followed  by  discussions  of  the  accuracy  of  the  pressure  solutions 
and  the  problem  of  inconsistency. 

Pressure  solutions  are  presented  as  contour  plots  of  lines 
of  constant  pressure,  or  isobars.  Referring  to  Fig.  4.12a,  the 
general  format  of  these  plots  is  as  follows: 

a.  The  dark  arrows  indicate  the  flow  direction. 

b.  All  solid  lines,  except  those  denoting  the  amplifier 
boundary,  are  lines  of  constant  static  pressure. 

c.  The  pressure  magnitudes  shown  relate  to  a reference 


value  of  one,  which  is  the  value  assigned  to  the  center 
node  of  the  inlet.  Due  to  the  non-dimensional ization , 
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a one  unit  change  in  magnitude  between  isobars  corre- 

2 

sponds  to  a dimensional  pressure  change  of  pll  . Water 

2 

flowing  at  73°F  at  a Reynolds  number  of  100  has  a pll  of 
2.2  x 10-4  in.  of  water,  or  8.01  x 10"^  lb/in^.  For  SAE  10 
oil  flowing  under  the  same  conditions,  this  quantity  is 
1.22  in.  of  water  or  4.39x10”^  lb^/in^. 

d.  The  dashed  lines  in  the  cavity  region  indicate  extensions 
of  the  isobars  into  regions  for  which  pressure  values 
were  not  calculated. 

Figures  4.12a  through  4.13b  are  pressure  fields  for  unde- 
flected flows  with  Reynolds  number  of  100.  Each  figure  differs  in 
the  ratio,  Q-j/Qs,  of  the  left  receiver  flow,  Q1 , to  the  supply  flow, 

Q . Figure  4.12a  represents  the  pressure  field  for  a symmetric  flow 
field.  The  pressure  drops  slightly  as  it  enters  the  cavity  region 
and  then  begins  a steady  increase  across  the  cavity.  The  pressure 
magnitude,  however,  varies  very  little  in  the  cavity  region  until 
the  flow  nears  the  splitter.  As  shown,  in  the  splitter  region  the 
pressure  gradients  are  relatively  large,  with  the  pressure  building 
to  a maximum  at  the  splitter.  The  pressure  then  decreases  equally 
in  both  receivers  as  the  flow  passes  out  of  the  amplifier. 

Similar  comments  can  be  made  about  Figs  4.12b  through  4.13b. 
In  each  successive  figure  of  this  series  more  flow  is  channeled 
through  the  left  receiver.  As  shown  in  Fig.  4.8a  through  Fig.  4.8d, 
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this  causes  the  central  jet  to  shift  to  the  left.  The  effect  on  the 
pressure  solution  is  to  move  the  region  of  highest  pressure  progres- 
sively farther  to  the  left,  corresponding  to  the  location  of  the 
main  jet  impingement  on  the  receivers. 

Figures  4.14a  and  4.14b  are  pressure  fields  for  deflected 
flows  which  have  Re =100,  equal  outflows,  and  0's  of  5°  and  15°,  re- 
spectively. For  these  cases  the  inlet  jet  itself  is  forced  into  the 

left  half  of  the  cavity  region.  This  results  in  a high  pressure  re- 
gion, as  shown  in  the  figures.  Because  the  outflows  are  equal,  the 
isobars  tend  back  towards  symmetry  at  the  splitter.  For  0=5°  the 
pressure  does  become  relatively  symmetric,  indicating  that  the  effect 
of  the  initial  deflection  is  small  on  the  receiver  region.  In  fact, 
the  pressures  at  the  receiver  outlets  are  0.574  and  0.563.  For 
0 = 15°,  however,  the  deflection  has  a noticeable  effect  in  the 

splitter  region,  with  the  higher  pressure  region  shifted  to  the  left 

as  would  be  expected.  Here  the  outlet  pressures  are  0.594  and  0.539. 

D . Discussion  of  Accuracy  of  Pressure  Solutions 

The  general  trends  of  these  pressure  solutions  seem  reason- 
able. There  is  some  question,  however,  as  to  the  accuracy  of  the 
pressure  magnitudes  predicted,  especially  for  the  case  having  one 
receiver  completely  blocked.  Several  approaches  were  taken  to  eval- 


uate the  accuracy  and  trends  of  the  pressure  magnitudes  obtained. 

The  first  approach  was  to  establish  an  upper  limit  on  the 
value  of  pressure  in  the  amplifier.  Consider  a flow  entering  the 


amplifier  with  a parabolic  profile.  Assume  that  there  are  no  viscous 

effects  on  the  system,  that  there  is  no  work  done  on  the  system,  and 

that  no  heat  is  transferred  to  the  system.  Assuming  that  the  flow 

moves  directly  across  the  amplifier  region,  the  center  streamline  in 

the  inlet  flow  would  impinge  on  the  splitter  tip.  At  the  inlet  the 

total  pressure  on  this  streamline,  in  non-dimensional  variables,  is 
2 

p + h v . Because  the  velocity  is  zero  on  the  splitter,  under  the 
assumptions  made  the  static  pressure  at  the  splitter  would  become 
equal  to  the  inlet  total  pressure.  In  the  absence  of  work  or  heat 
transfer,  the  introduction  of  viscous  effects  can  do  nothing  but 
lower  this  value.  Therefore,  the  highest  total  pressure  in  the  inlet 
forms  an  upper  limit  on  the  magnitude  of  the  static  pressure  in  the 
amplifier.  For  this  amplifier  this  limit  would  be  p = 2.125. 

The  highest  pressure  in  the  receiver  inlet  region  for  the 
symmetric  case,  p = 1.616,  certainly  fits  within  this  limit.  In 
fact,  as  listed  in  Table  4.2,  the  highest  pressure  for  all  cases, 
except  the  blocked  case,  meets  this  criterion.  It  can  be  concluded 
then  that  the  magnitude  for  the  blocked  case  is  incorrect  but  that 
for  all  other  cases  the  values  are  at  least  possible. 

Further  examination  of  the  pressure  magnitudes  can  be  made 
by  using  the  line  integral  method  discussed  earlier.  This  method 
can  be  implemented  easily  for  a straight  path  between  the  center  of 
the  inlet  and  the  tip  of  tr.e  splitter.  Eq . (2.9)  was  numerically 
integrated  usi-ng  Simpson's  rule  [1]  with  the  necessary  derivatives 
supplied  through  centered  difference  approximations  involving  the 
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Pressure  on  Splitter 


Ql/Qs 

Highest 

Pressure 

Poisson 

Solution 

Line 

Integral 

Inconsistency 

0.000 

3.647 

3.340 

1.235 

2.039 

0.250 

2.119 

2.093 

1.367 

0.446 

0.375 

1.750 

1.750 

1.348 

0.002 

0.500 

1.616 

1.616 

1.388 

-0.171 

0.625 

1.750 

1.750 

1.348 

0.002 

0.750 

2.119 

2.093 

1.367 

0.446 

1.000 

3.647 

3.340 

1.235 

2.039 

Table  4.2.  Pressure  Solution  Details. 

All  quantities  non-dimensional. 

Re  = 100. 

flow  field  solutions.  For  the  symmetric  case  the  value  obtained  for 
the  pressure  at  the  splitter  tip  is  1.388.  This  is  15*  lower  than 
the  value  predicted  by  the  solution  of  Eq.(2.8),  which  will  be  re- 
ferred to  as  the  Poisson  solution.  As  discussed  earlier,  error  is 
possible  in  use  of  the  line  integral  method.  Which  answer  is  the 
better  is  still  in  question.  The  point  is  that  these  values  are  cer- 
tainly comparable  and  that  either  value  is  reasonable.  It  should  be 
emphasized  again  that  both  methods  utilize  identical  flow  field  solu- 
tions and  any  differences  here  are  due  entirely  to  the  different 
methods  for  extracting  pressures. 

As  the  flow  rate  in  the  receivers  is  varied  from  the  symmet- 
ric case,  the  line  integral  method  predicts  very  little  change  in 
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pressure,  as  tabulated  in  Table  4.2.  In  comparison  the  values  pre- 
dicted by  the  Poisson  solution  become  sharply  larger. 

Some  indication  of  which  trend  is  correct  can  be  obtained  by 
a control  volume  analysis  of  the  central  jet,  similar  to  that  which 
forms  the  basis  for  the  cowl  streamline  method  [11].  The  results  will 
only  illustrate  trends  because  of  the  necessary  simplifying  assump- 
tions. These  are 

a.  The  viscous  loss  across  the  cavity  is  negligible. 

b.  The  pressure  across  both  the  inlet  and  the  splitter 
regions  is  constant. 

It  is  also  noted  that  there  is  no  work  done  or  heat  transferred  to 
the  flow.  The  control  volume  is  taken  as  the  region  between  the 
streamlines  t|i  = 0 and  ip  = -1.0.  For  the  symmetric,  Re  = 100  case 
this  would  result  in  the  control  volume  shown  in  Fig.  4.15.  For  the 
assumed  flow  the  applicable  non-dimensional  energy  equation  [27]  is 


I (H  V2  + p)  vdx 
inlet 


ih  + p) 
receiver 


(4.1) 


? 2 

where  V = / v + u . Using  the  velocities  given  by  the  numerical 

flow  solutions  and  assuming  the  pressure  at  the  inlet  equal  to  one, 
Eq. (4.1)  can  be  solved  for  the  average  pressure  over  the  receiver  end 
of  the  control  volume.  The  integrals  were  evaluated  using  Simpson's 
rule. 


For  the  symmetric  case  the  end  of  the  control  volume  extends 
across  both  receivers.  The  average  pressure  predicted  is  0.785.  For 
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Figure  4.15.  Typical  Control  Volume  for 

Integral  Pressure  Calculations 

the  blocked  case  the  control  volume  covers  only  the  left  receiver. 

For  this  region  the  predicted  average  pressure  is  1.305.  Because 
these  pressures  are  averages,  no  comparisons  will  be  made  using  the 
actual  magnitudes  predicted.  The  important  comparison  is  that  there 
is  a substantial  increase  in  the  average  pressure  in  the  receiver 
inlet  region  as  the  outflow  increases  through  the  left  receiver. 

This  general  trend  is  the  same  as  that  found  in  the  Poisson  solutions 
for  pressure  and  is  contrary  to  that  predicted  by  the  line  integral 
method.  In  this  respect  then  the  Poisson  solutions  are  preferred. 

A further  check  on  the  accuracy  of  the  pressure  solutions 
can  be  obtained  by  examination  of  the  predicted  pressure  drop  in  the 
receivers.  At  the  outlet  of  the  two  receivers  the  flow  has  been 
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forced  to  be  fully  developed  by  the  imposed  boundary  conditions. 
Moreover,  with  this  restriction  at  the  boundary  the  resulting  flow 
solutions  predict  that  the  flow  is  fully  developed  over  the  entire 
last  half  of  each  receiver.  Hence,  this  flow  is  fully  developed 
flow  between  parallel  plates,  for  which  there  is  an  analytic  solu- 
tion for  the  pressure  drop.  In  terms  of  non-dimensional  variables 
the  pressure  gradient  is  [27] 

dP  = _ 1?  (Ml  (4.2) 

da  wJ  Re 

where  Aip  is  the  difference  between  \p  on  the  two  walls.  Evaluation 
of  Eq. (4.2)  gives  an  analytic  prediction  of  the  pressure  gradient 
for  a given  flow  situation.  A,  numerical  prediction  can  be  obtained 
by  centered  differencing  of  the  predicted  pressure  values  in  this 
region.  For  the  symmetric,  Re  = 100  case,  Eq. (4.2)  predicts  a pres- 
sure gradient  of  -0.1697  for  the  left  receiver.  Differentiation  of 
the  numerical  results  gives  -0.1654.  Similarly  for  the  cases  of 
Q-j  = 0.75  Qs  and  Q1  = 1.0  Q , the  comparisons  are,  respectively, 
-0.254  versus  -0.264  and  -0.340  versus  -0.386. 

The  most  obvious  disagreement  occurs  in  a blocked  receiver. 
Here  there  should  be  no  pressure  gradient  at  all,  but  a gradient  of 
-0.046  is  predicted  numerically.  This  is  caused  by  the  method  for 
compensation  for  inconsistency.  The  compensation  introduces  a non- 
zero s..  into  the  equations  for  this  region.  Nevertheless,  the 
• J 

general  agreement  is  good  and  this  comparison  indicates  that  the 
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Poisson  solution  does  indeed  yield  the  correct  pressure  gradients  in 
the  receivers  for  the  given  flow  fields. 

Because  experimental  pressure  values  were  not  obtained,  the 
above  comparisons  form  the  basis  for  an  evaluation  of  the  pressure 
results.  As  stated,  the  location  of  the  high  pressure  regions  seems 
to  be  correct  for  the  different  cases.  Also,  for  the  symmetric  case 
the  pressure  magnitudes  are  reasonable.  As  the  flow  conditions  vary 
away  from  symmetric  flow,  the  pressure  values  in  the  region  of  high- 
est pressure  change  in  the  correct  manner,  i.e.  they  increase  over 
the  more  symmetric  case.  However,  they  increase  too  much  and  for 
the  blocked  case  reach  magnitudes  which  are  not  physically  possible. 

In  summary,  it  appears  that  results  of  the  numerical  model  are  reason- 
able for  cases  in  the  approximate  range  of  0.250  Q$  < Q-|  < 0.750  Q$. 
Outside  of  this  range  the  results  are  questionable. 

A possible  explanation  for  this  behavior  involves  the  wall 
gradients  and  internal  derivatives  of  which  are  being  numerically 
evaluated.  For  the  symmetric  case  half  of  the  flow  leaves  through 
each  receiver.  As  the  specified  flow  situation  deviates  from  the 
symmetric  case,  more  flow  must  exit  through  one  receiver  until,  for 
the  blocked  case,  all  the  flow  exits  through  one  receiver.  As  a 
result  of  increased  flow  through  the  same  area,  the  gradients  in 
stream  function  in  the  region  near  the  splitter  must  necessarily 
become  larger.  Further,  these  gradients  must  be  evaluated  for  use 
in  the  Poisson  solution  using  finite  difference  approximations  hav- 
ing the  same  grid  spacing.  Apparently,  as  the  flow  rate  increases 
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the  five  grid  points  across  each  receiver  are  not  sufficient  to 
allow  accurate  determination  of  these  gradients.  When  these  in- 
accurate gradients  are  used  in  the  Poisson  equation,  the  previously 
discussed  inaccuracies  in  the  pressure  values  result. 

Besides  the  question  of  inaccuracies  discussed  above,  there 
remains  the  problem  of  inconsistency.  A physical  interpretation  of 
inconsistency  can  be  obtained  by  examining  a two  dimensional  heat 
conduction  problem.  The  appropriate  equation  is 


(4.3) 


where  T is  the  temperature,  u"'  is  the  internal  heat  generation,  and 
k is  the  thermal  conductivity  [2].  Consider  a square  region  with 
insulation  on  all  four  sides.  The  heat  transfer  across  the  bound- 
aries is  zero  and  the  resulting  boundary  conditions  are 


9T 

3n 


0 


(4.4) 


For  this  Neumann  problem  consistency  requires  that  the  integral  of  the 
forcing  function,  um/k,  over  the  region  must  equal  the  line  integral 
of  the  normal  temperature  gradient  on  the  boundary,  that  is, 


dxdy  = <p 


( 3T 


3n 


ds 


(4.5) 


Physically  this  means  that  whatever  heat  is  generated  in  the  region 
must  either  be  absorbed  in  the  region  or  must  leave  the  region. 
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Because  no  heat  can  leave  the  region  for  the  problem  considered, 
the  integral  of  the  internal  heat  generation  must  be  zero. 

Other  heat  conduction  problems  can  be  imagined  to  further 
illustrate  the  concept  of  inconsistency.  These  problems  are  easily 
visualized  in  physical  terms  and,  consequently,  they  form  a conven- 
ient method  for  viewing  inconsistency. 

The  physical  interpretation  is  not  nearly  as  simple  for  the 

pressure  problem.  It  is  not  clear  what  is  physically  represented  by 

'2  2 2,2] 

the  forcing  function  of  Eq.(2.8),  which  is  2 . How- 

lax2  9y2  3x^yJ 

ever,  this  function  is  mathematically  analogous  to  the  heat  genera- 
tion term  in  the  conduction  problem.  Therefore,  the  boundary 
gradients  on  pressure  must  allow  for  this  "generation. " If  they  do 
not,  the  problem  is  inconsistent  and  is  a contradiction  of  physical 
principles. 

In  contrast,  for  the  corresponding  difference  equations  a 
converged  solution  can  be  obtained  through  iteration  using  SOR.  This 
assumes  the  convergence  criterion  used  in  this  study.  As  pointed  out 
in  Section  II-C,  the  numerical  algorithm  appears  to  allow  for  the 
inconsistency  as  it  converges.  The  mathematical  reason  and  explana- 
tion for  this  was  not  found  but  a series  of  numerical  experiments 
were  performed  in  an  effort  to  characterize  the  phenomenon. 

Besides  the  previously  discussed  residuals,  the  most  interest- 
ing effect  observed  dealt  with  the  convergence  of  an  inconsistent 
formulation.  Using  a conventional  convergence  criterion,  say 
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k+1  k 

- P ^ j | < e , with  no  manipulation  of  the  pressure  arrays,  as 
in  Section  II-C,  the  solution  may  or  may  not  converge.  The  differ- 
ences between  successive  iterations  diminish  until  the  pressure  at 

k+1  k 

each  point  varies  a fixed  amount,  that  is,  i P - P - = constant. 
The  magnitude  of  this  change  depends  on  the  value  of  the  inconsis- 
tency, the  relaxation  factor,  and  the  grid  spacing.  If  the  inconsis- 
tency is  large  enough,  the  change  will  never  become  less  than  . 

Even  if  the  convergence  criterion  is  met,  the  solution  values  wou^ 
still  drift  if  the  calculations  were  continued.  Thus,  consideration 
of  inconsistency  helps  explain  the  drifting  behavior  noted  cy  Roache 
[20]  for  convergence  of  pressure  solutions.  Further,  it  was  observed 

that  if  the  pressure  formulation  were  made  consistent,  the  solution 

I k+1  k ' 

would  indeed  converge  and  the  change  P.j  - P^j  would  go  to  zero. 

The  convergence  criterion  used  in  this  study  is  not  affected 
by  the  above  problem  because  it  forces  convergence  to  a surface  shape 
and  is,  therefore,  not  affected  by  uniform  translation  of  the  entire 
surface.  It  is  felt  that  such  a criterion  is  the  appropriate  con- 
vergence criterion  for  the  Neumann  problem. 

The  effect  of  inconsistency  on  different  iterative  solution 
techniques  was  not  determined.  All  observations  were  made  using 
point  SOR  as  the  method  of  solution.  One  other  method,  Jacobi's 
method  [20,30],  was  attempted  on  a simple  Neumann  problem.  The 
resulting  solution  exhibited  an  apparent  standing  wave  propagating 
back  and  forth  across  the  region.  Based  on  this  experience. 
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Jacobi's  method  seems  inappropriate  for  use  on  the  Neumann  prob- 
lem. 

2 

Convergence  of  V P for  an  inconsistent  formulation  to  the 
values  s..  + s'  naturally  results  in  an  incorrect  relationship 
between  values  of  pressure  at  adjoining  nodes.  In  particular,  the 
gradients  specified  as  boundary  conditions  are  not  reflected  in 
the  converged  solutions.  This  observation  is  also  made  by  Young 
and  Gregory  [30]. 

The  actual  magnitude  of  inconsistency  varies  between  the 
different  solutions.  For  Re  = 100  the  symmetric  case  had  a total 
inconsistency  of  -0.171.  The  effect  of  this  inconsistency  is  shown 
rather  dramatically  by  comparing  Fig.  4.16a  and  Fig.  4.16b.  Here 
isobars  of  equal  magnitude  are  plotted  for  compensated.  Fig.  4.16a, 
and  non-compensated.  Fig.  4.16b,  pressure  solutions  based  on  the 
same  flow  field.  It  should  be  noted  that  this  comparison  is  some- 
what misleading,  since  the  large  change  in  the  isobars  was  caused 
by  a small  change  in  the  pressure  values.  An  indication  of  the 
magnitude  of  change  between  the  respective  solutions  is  given  by 
comparison  of  the  value  predicted  for  pressure  at  the  splitter  tip. 
These  values  are  1.616  for  the  compensated  case  and  1.575  for  the 
non-compensated  case. 

Referring  to  Table  4.2,  the  magnitude  of  inconsistency 
becomes  smaller  for  the  case  where  = 0.625  Q$  and  then  increases 
as  the  blocked  case  is  approached.  This  follows  a similar  trend  as 
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Figure  4.16.  Comparison  of  Compensated  and 

Non-compensated  Pressure  Solutions 
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the  inaccuracy  in  pressure  magnitudes  discussed  earlier  and  the  two 
problems  are  probably  related  through  the  inaccurate  resolution  of 
gradients  near  the  splitter. 

It  is  felt  that  there  are  two  probable  causes  of  numerica' 
inconsistency.  Both  are  related  to  truncation  error.  The  first,  and 
most  obvious,  is  the  inaccuracies  inherent  in  forming  the  difference 
approximations  for  the  boundary  gradients  and  the  forcing  functions. 
This  is  the  same  problem  as  discussed  earlier  as  an  explanation  for 
the  inaccuracy  of  the  pressure  magnitudes.  Given  an  accurate  flow 
solution  and  a means  for  accurate  evaluation  of  the  derivatives  of  > p 
in  Eqs . (2.8)  and  (2.9),  there  would  be  no  inconsistency  problem.  The 
flow  equations  are  derived  taking  the  pressure  gradients  into  account. 
Therefore,  if  flow  solutions  were  exact,  the  gradients  and  forcing 
functions  would  necessarily  be  compatible.  However,  this  is  not  the 
case.  There  is  some  error  in  the  flow  solutions  and,  more  signifi- 
cantly, there  is  the  possibility  of  relatively  large  errors  in  the 
necessary  numerical  derivatives  of  Again,  this  is  compounded 
when  the  gradients  become  steeper  and  more  difficult  to  approximate 
as  Q-|  increases.  Such  errors  in  the  approximations  result  in  a fail- 
ure to  satisfy  the  discrete  consistency  equation. 

Although  there  is  some  error  in  derivative  approximation 
throughout  the  region  it  is  felt  that  the  bulk  of  the  inconsistency 
results  from  errors  in  particular  regions.  If  the  difference  approxi- 
mations could  be  improved  in  these  areas,  the  inconsistency  problem 
should  be  greatly  alleviated. 
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A second  possible  cause  of  numerical  inconsistency  lies  in  the 
conservation  properties  of  the  basic  finite  difference  scheme  for  the 
flow  field.  Second  upwind  differencing  is  designed  to  conserve  ver- 
ticity  [23].  By  this  it  is  meant  that  in  any  region  the  time  rate 
of  accumulation  of  s equals  the  net  advection  and  Diffusion  jf  £ 
across  the  boundaries  of  the  region  [20].  This  is  a property  of  the 
difference  scheme  and  can  be  guaranteed  by  proper  formulation  of  the 
difference  approximations  to  Eq. (2.6) . Other  schemes  are  designed  to 
conserve  the  vorticity  squared,  linear  momentum,  and/or  kinetic 
energy  [20].  It  is  speculated  that  there  is  a quantity  related  to 
pressure,  which  may  or  may  not  be  one  of  the  above,  which  also  must 
be  conserved  by  the  differencing  scheme.  Possibly  this  quantity 
could  be  isolated  and  differencing  schemes  found  in  the  literature, 
or  developed,  which  would  give  the  desired  conservation  properties. 

E . Characteristic  Design  Curves 

Recall  that  for  design  purposes  a steady  flow  fluidic  ampli- 
fier is  described  by  one  or  more  characteristic  curves.  These  can 
be  computed  from  the  mathematical  model  by  selecting  ar.  inlet  flow 
rate,  or  Reynolds  number,  and  a deflection  <j>.  The  flow  rates  in  the 
receivers  can  then  be  varied  from  blocked  to  full  flow  and  a flow 
field  and  resulting  pressure  solution  obtained  for  each  case.  T).e 
pertinent  pressure  values  at  the  amplifier  ports  can  be  plotted  on 
any  of  the  various  characteristic  curves  described  in  Chapter  I. 
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For  example,  the  flow  in  one  receiver  can  be  plotted  against  the 
pressures  at  the  outlet  of  the  receivers,  as  in  Fig.  4.17. 

One  series  of  solutions  for  positive  <p  will  also  give  the 
results  needed  for  negative  <p  due  to  symmetry  conditions.  That  is, 

P1  for  the  case  <j>  = 15°,  and  Q-j  = 0.625  Qs  is  the  same  as  P^  for 
(J)  = -15°,  Q2  = 0.625  Q$.  Similarly,  P-j  for  the  case  <J>  = 0°  and 
Q-j  = 0.625  Q$  is  the  same  as  P^  for  the  case  <J>  = 0°  and  Q2  = 0.375 
Qs<  Thus,  only  eight  solutions  were  required  to  generate  the  15 
points  used  in  Fig.  4.17. 

The  accuracy  of  the  curves  plotted  in  Fig.  4.17  for  the  device 
in  this  study  are  questionable  due  to  the  problems  in  obtaining  accu- 
rate pressure  solutions.  These  were  discussed  in  the  previous  section. 
Because  of  the  accuracy  problem  only  the  cases  with  receiver  flow 
rates  between  0.250  Q$  and  0.750  Q$  were  plotted.  Even  with  this  re- 
striction, the  trends  do  not  seem  quite  right  compared  with  physical 
reasoning. 

For  example,  the  double  valued  character  around  P = 0.5  is 
unlikely.  Consider  the  curves  for  P^.  This  outlet  pressure  is 
determined  by  the  pressure  in  the  receiver  inlet,  the  pressure  drop 
due  to  entrance  effects  as  the  flow  enters  the  receiver,  and  the 
viscous  pressure  drop  in  the  receiver.  The  pressure  in  the  splitter 
region  is  expected  to  rise  as  the  flow  increases  (Section  D,  Chap- 
ter IV).  Exactly  what  the  entrance  effects  would  do  as  the  flow  rate 
increases  is  uncertain.  The  steady  state  viscous  pressure  drop  in 
the  receivers  increases  linearly,  as  in  Eq.(4.2).  Using  the  pressure 
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gradients  presented  in  Section  D and  a length  of  5 units,  the  linear 
pressure  drop  in  the  receiver  would  increase  from  roughly  -0.827  for 
the  symmetric  case  to  -1.93  for  the  blocked  case.  For  this  same 
range  of  flow  rates  the  receiver  inlet  pressure  begins  around  1.6 
and  is  bounded  by  an  upper  value  of  2.125.  Since  entrance  effects 
must  increase  the  pressure  drop  down  the  receiver,  it  is  felt  that 
the  terms  producing  a pressure  drop  in  the  receiver  should  dominate 
the  pressure  increase  at  the  inlet  and  the  curve  for  should  con- 
tinue to  lower  values  of  pressure  as  0^  increases. 

It  is  also  noted  that  the  cross-over  of  the  curves  for  the 
deflected  cases  at  approximately  P = 0.75  is  questionable.  Neverthe- 
less, with  refinement  of  the  pressure  solutions  there  is  no  reason 
why  useful  characteristic  curves  for  design  work  could  not  be  ob- 
tained for  the  oeometry  considered. 

F.  Edgetone  Oscillations 

For  a laminar  free  jet  impinging  on  a solid  object  it  has 
been  noted  in  the  literature  [4,5,22]  and  was  observed  in  the  experi- 
ments of  this  study  that  a transverse  oscillation  of  the  jet  occurs 
for  certain  flow  situations.  These  "edgetone"  or  "jet-edge"  oscilla- 
tions are  stable,  periodic  and  two-dimensional. 

An  example  of  this  oscillation  observed  in  the  present  study 
is  shown  in  Fig.  4.18.  For  this  flow  the  left  receiver  is  blocked 
and  the  Reynolds  number  is  360.  Each  picture  in  Fig.  4.18  represents 
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the  same  streamline  at  a different  point  in  time.  The  oscillation 
was  stable,  had  a constant  frequency,  f,  of  0.22  cycles/sec,  and 
was  of  relatively  large  amplitude. 

The  first  signs  of  oscillation  were  experimentally  observed 
before  the  Reynolds  number  reached  360.  At  Re = 290  the  first  signs 
of  an  oscillation  were  observed  and  when  Re = 300  the  oscillation  was 
well  defined.  This  oscillation  had  a slightly  lower  frequency, 

0.19  cycle  per  second,  and  a much  lower  amplitude  than  the  oscilla- 
tion shown  in  Fig.  4.18. 

As  discussed  by  Rockwell  [22],  edgetone  oscillations  depend 
on  the  Reynolds  number,  the  distance  h the  incoming  jet  must  travel 
before  impinging  on  a surface,  the  width  b of  the  inlet  nozzle,  and 
the  height  of  the  cavity.  The  effect  of  cavity  height  is  to  intro- 
duce three  dimensional  effects  which  tend  to  retard,  with  respect 
to  increasing  Reynolds  number,  the  onset  of  oscillation  [22,25]. 

In  his  experimental  work,  which  was  highly  three  dimensional  since 
the  nozzle  aspect  ratio  equaled  one,  Rockwell  observed  that  no  oscil- 
lations occurred  in  flows  which  had  Reynolds  numbers  less  than  1250 
and  a h/b  ratio  of  10.  With  the  nozzie  aspect  ratio  of  4 used  in 
the  experiments  of  this  study  the  three  dimensional  effects  are 
reduced  and  the  Reynolds  number  at  the  onset  of  oscillation  should 
be  reduced,  as  was  the  case. 

In  the  present  study  the  value  for  the  Strouhal  number,  fb/U, 
was  identical  for  the  two  experimental  oscillations  for  which  f was 
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measured.  This  value  was  0.0095.  That  it  is  a constant  agrees 
with  the  trends  presented  by  Rockwell  [22].  The  magnitude,  however, 
is  lower  than  that  predicted  by  Rockwell.  An  extrapolation  of 
Rockwell's  results  to  the  flow  conditions  of  this  study  predict  a 
Strouhal  number  value  of  0.03. 

Unfortunately,  the  edgetone  oscillations  were  not  predicted 
by  the  numerical  simulation.  Therefore,  although  the  simulation  is 
stable  and  will  converge  for  higher  Reynolds  numbers,  the  results  are 
meaningful  only  up  to  a Reynolds  number  of  roughly  290.  A probable 
cause  of  the  failure  of  the  numerical  model  to  predict  the  edgetone 
oscillations  is  the  handling  of  the  outflow  and  inflow  boundary  con- 
ditions. In  both  cases  the  flow  is  forced  to  be  fully  developed. 

In  effect,  an  oscillation  cannot  exist  at  the  receiver  outlet  nor 
at  the  amplifier  inlet.  These  constraints  tend  to  suppress  the 
development  of  oscillations  in  the  flow  field. 

A second  possible  cause  of  this  weakness  in  the  numerical 
model  is  the  artificial  viscosity  introduced  by  second  upwind  differ- 
encing [20,21].  Artificial  viscosity  refers  to  a completely  numeri- 
cal, artificial  diffusion  effect  which  is  introduced  by  the  method 
of  differencing.  This  term  damps  out  small  perturbations  and  in- 
stabilities in  the  flow  and  could  be  eliminating  the  physical 
mechanism  which  triggers  the  edgetone  oscillations. 
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G . Discussion  o f the  Numerical  Met hod 

During  the  development  of  the  numerical  model  for  the  flow 
field  several  observations  on  the  computational  procedure  itself  were 
made  beyond  those  mentioned  earlier  in  this  chapter.  These  observa- 
tions will  be  presented  in  this  section. 

Discussion  in  Section  D of  this  chapter  emphasized  the 
possible  resolution  problems  in  the  splitter  region.  The  possibility 
of  resolution  problems  was  considered  even  before  these  problems  were 
discovered  in  the  pressure  solution  work.  Early  in  the  study  compari- 
sons were  made  between  solutions  to  the  same  case  but  with  two  sep- 
arate grids.  One  grid  consisted  of  21  nodes  in  the  x-direction  and 
29  nodes  in  the  y-direction,  giving  Ax  = Ay  = 0.5.  The  other  grid 
was  the  41  x 57  grid  used  for  all  results  presented  in  this  report. 

It  has  Ax  = Ay  = 0.25. 

For  the  lower  Reynolds  numbers  the  results  on  the  different 
grids  were  very  similar.  For  example,  compare  Fig.  4.19a  and  Fig. 
4.19b,  which  are  such  a set  of  solutions  for  symmetric  flow  with 
Re =10.  However,  for  larger  values  of  Reynolds  number  the  results 
were  different  on  the  two  grids.  Figure  4.19c  and  Fig.  4 . 1 9d  are 
examples  of  such  a difference  for  symmetric.  Re =100  flow. 

This  result  is  reasonable  when  the  accuracy  of  the  difference 
scheme  is  considered.  Recall  that  the  differencing  of  the  vorticity 
equation  is  second-order  accurate  on  the  differencing  of  the  diffu- 
sion terms  but  only  first-order  accurate  on  the  differencing  of  the 


21  x 29  Grid  d)  Re  = 100  41  x 57  Grid 

Figure  4.19.  Effect  of  Grid  Spacing 
on  Flow  Solutions 
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advection  terms.  For  low  Reynolds  numbers  the  flow  is  diffusion 
dominated  and  the  second-order  accurate  differencing  gives  similar 
solutions  for  both  grids.  At  higher  Reynolds  numbers  the  advection 
terms  dominate  the  flow  field  and  the  first-order  accurate  differ- 
encing is  insufficient  to  model  the  process  on  the  coarse  grid. 

This  results  in  the  large  difference  between  the  two  solutions  for 
Re  = 100. 

As  discussed  earlier  the  experimental  flow  fields  and  those 
predicted  using  a 41  x 57  grid  compare  very  well.  Conversely,  the 
pressure  results  indicate  that  regardless  of  the  agreement  between 
experimental  and  numerical  flow  results  the  grid  needs  refining  in 
the  splitter  region  to  obtain  good  pressure  solutions.  Thus,  further 
resolution  studies  should  be  aimed  at  examining  such  techniques  as 
variable  grid  spacing  so  that  extra  nodes  are  placed  near  the 
splitter. 

As  is  probably  true  in  most  computational  fluid  dynamics 
problems,  the  final  numerical  boundary  conditions  were  not  the  same 
t those  initially  attempted.  Two  changes  worth  noting  are  illus- 
trated in  Fig.  4.20.  All  three  contour  plots  represent  converged 
flow  fields  for  symmetric  flow  with  Re =100.  Figure  4.20a  is  the 
first  converged,  stable  solution  for  Re=100  obtained  in  this  study. 
In  this  formulation  the  wall  vorticity  in  the  cavity  region  was 
given  by 
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where  n is  the  normal  to  the  wall  [20].  Note  that  this  expression  is 
first  order  accurate.  This  equation  was  also  used  to  evaluate  the 
vorticity  at  the  corners  of  the  inlet  instead  of  setting  these  two 
values  in  keeping  with  the  parabolic  inlet  velocity  profile  (Chapter 
II,  Section  C).  Because  the  vorticity  is  singular  at  these  two 
points,  it  was  not  initially  apparent  which  of  these  alternatives 
to  take. 

The  first  refinement  on  the  wall  vorticity  resulted  in  Fig. 
4.20b.  The  formulation  differs  in  that  the  second  order  accurate 
expression  for  wall  vorticity,  Eq.(2.18e),  was  implemented.  Notice- 
able change  occurred  in  the  recirculation  regions,  with  the  maximum 
value  of  >p  in  the  main  recirculation  region  becoming  smaller  and 
the  secondary  recirculation  region  disappearing  altogether.  However, 
the  general  shape  of  the  flow  remained  the  same  and,  in  fact,  the 
main  flow  was  not  changed  at  all. 

The  final  change,  resulting  in  the  flow  field  of  Fiq.  4.20c, 
was  to  try  the  other  alternative  for  the  inlet  corner  vorticity 
values.  Instead  of  using  Eq . (2.1 8e ) these  values  were  based  on  the 
parabolic  inlet  velocity  profile,  as  in  Eq.(2.20).  The  effect  on 
the  main  jet  is  clearly  visible.  The  0.0  and  -1.0  streamlines  detach 
from  the  wall  immediately  and  extend  downstream.  The  distance  be- 
tween the  two  streamlines  remains  considerably  less  than  in  (b). 

That  this  formulation  was  the  better  choice  for  a boundary  condition 
is  confirmed  by  the  comparisons  with  experimental  results. 
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An  important  consideration  in  the  use  of  a numerical  routine 
as  a design  tool  is  the  computational  time  required  to  obtain  a con- 
verged solution  as  this  translates  directly  into  costs.  For  several 
reasons  it  is  difficult  to  give  a meaningful  figure  for  convergence 
time.  In  particular,  the  same  numerical  routines  would  require 
varying  times  for  convergence  on  different  computer  systems.  The 
figures  given  here  will  all  relate  to  the  CDC  6600  system  in  use  at 
The  University  of  Texas.  Even  on  that  particular  system  several 
numbers  could  be  used  as  the  figure  for  convergence  time.  The  time 
presented  here  is  the  "TM"  time,  which  includes  central  processing 
time,  peripheral  processing  time,  and  a time  charge  based  on  tape 
and  disk  usage. 

As  an  example  consider  the  symmetric  case  for  flow  with 
Re =100.  The  flow  field  for  this  case  converged  in  756  seconds. 

The  pressure  solution  then  required  28  seconds. 

The  time  required  for  the  pressure  solution  varied  very  little 
for  the  different  cases.  The  time  for  obtaining  the  flow  field  solu- 
tion varied  considerably  for  the  different  cases,  with  the  variations 
following  some  general  trends.  As  the  Reynolds  number  increased  the 
time  required  for  convergence  increased.  For  symmetric  flow  with 
Reynolds  numbers  of  20,  100,  and  275,  these  times  were  respectively 
693,  756,  and  1021  seconds.  Although  the  exact  time  for  Re =1000 
is  not  available,  that  solution  required  roughly  2500  seconds. 

Also,  for  a given  Reynolds  number  unsymmetric  flows  in 
general  required  more  time  for  convergence.  Consider  the  effect 


of  unsymmetric  outflow.  For  Re =100  and  no  deflection,  the  symmetri 
case  required  756  seconds  while  the  blocked  case  required  1528  sec- 
onds. Likewise,  with  balanced  outflows  but  with  a 15°  deflection 
of  the  inlet  jet,  the  convergence  time  was  1058  seconds.  These 
trends  seem  reasonable  because  the  flow  becomes  physically  more  com- 
plex as  the  case  considered  becomes  more  unsymmetric.  However,  it 
should  be  noted  that  there  were  deviations  from  these  trends. 

Consider  the  computer  time  required  to  produce  the  character 
istic  curves  shown  on  Fig.  4.17.  The  fifteen  points  presented  on 
this  figure  required  solutions  to  eight  separate  flow  cases.  The 
flow  field  and  pressure  solutions  for  these  eight  cases  took  a total 
of  6768  seconds  to  compute,  which  is  an  average  of  846  seconds  per 
case. 

At  this  stage  in  the  development  of  the  computer  model  as  a 
design  tool  this  figure  indicates  that  such  a use  might  be  competi- 
tive with  other  approaches.  However,  no  actual  analysis  was  made  to 
determine  what  an  "economic"  cost  for  a simulation  might  be.  It  is 
simply  noted  that  the  actual  computer  time  is  not  unduly  high, 
especially  considering  the  following  points. 

First,  the  computer  time  is  directly  related  to  the  degree 
of  convergence  required  of  the  routines,  and  the  solutions  obtained 
have  a high  degree  of  convergence.  A short  study  was  performed  to 
determine  the  effect  of  the  convergence  of  the  flow  field  on  the 
pressure  solution.  For  the  symmetric,  Re =100  case  pressure 


solutions  were  obtained  at  intermediate  time  steps  109,  144,  and  181 
in  addition  to  the  final  "converged"  solution  at  time  step  204. 

These  pressure  solutions  were  compared  to  see  what  was  the  absolute 
change  in  pressure  between  the  four  times.  Between  time  steps  109 
and  144  the  pressure  changed  significantly,  with  a maximum  change  of 
0.0837  among  pressure  values  with  a magnitude  around  1.0.  Between 
steps  144  and  181  and  between  steps  181  and  204,  however,  the  maxi- 
mum change  in  pressure  dropped  to  0.0108  and  0.0018  respectively. 
Thus,  it  could  conservatively  be  said  that  the  pressure  was  within 
2%  of  its  final  value  after  time  step  144  in  the  flow  field  solution. 
This  corresponds  to  a relative  convergence  criterion  on  vorticity  of 
0.008,  rather  than  the  0.001  used  in  this  study.  It  also  corresponds 
to  a 23%  savings  in  computation  time.  Whether  such  a savings  could 
be  realized  on  all  cases  is  not  known,  but  it  is  strongly  felt  that 
a less  stringent  error  criterion  could  be  used  without  sacrificing 
significant  quality  in  the  solutions. 

Second,  beyond  the  question  of  a convergence  criterion  there 
is  room  for  improvement  in  the  convergence  times  in  both  the  numeri- 
cal method  used  and  the  coding  of  the  method.  Although  the  program 
was  written  with  speed  in  mind,  there  is  certainly  possible  improve- 
ment in  this  respect  since  it  was  not  written  by  a programming 
specialist.  Further,  the  use  of  quicker  but  more  complex  solution 
methods  could  lower  convergence  times.  Either  or  both  of  these 
improvements  could  be  made  if  the  convergence  speed  became  a critical 
factor  in  the  use  of  the  model  as  a design  tool. 


The  final  observation  of  this  section  deals  with  the  tran- 
sient behavior  of  the  numerical  solution.  Recall  that  the  time 
dependent  vorticity  transport  equation  was  used  as  the  basic  equation 
for  the  flow  field.  Also  recall  that  the  n+1  boundary  values  on  wall 
vorticity  are  calculated  and  used  at  each  time  step.  Thus,  given 
that  the  time  step  is  kept  small  enough  to  produce  accurate  solu- 
tions, there  is  no  reason  that  the  time  dependent  behavior  of  the 
flow  field  cannot  be  predicted  as  the  flow  field  changes  from  some 
initial  state  to  some  later  state.  In  fact,  by  examining  several  of 
the  intermediate  flow  solutions  for  a symmetric.  Re =1000  case  it  was 
observed  that  disturbances  were  convected  through  the  amplifier  in  a 
time  like  manner. 

There  are,  however,  several  reservations  to  acceptance  of 
this  as  true  transient  flow  First,  the  inability  of  the  model  to 
predict  the  edgetone  oscillations  indicates  that  some  unsteady,  or 
transient,  flow  phenomena  are  not  being  modeled,  as  discussed  in 
Section  F.  Another  unresolved  problem  is  what  should  be  the  physi- 
cal interpretation  of  the  initial  conditions.  For  example,  does 
using  the  symmetric.  Re =100  converged  solution  as  the  initial  flow 
distribution  and  changing  the  Reynolds  number  from  100  to  200  corre- 
spond to  a step  increase  in  the  average  inflow  velocity,  U?  Although 
a tempting  interpretation,  such  an  interpretation  is  not  completely 
consistent  w:th  the  assumptions  made  in  deriving  the  basic  equa- 
tions. In  particular,  the  non-dimensional izing  assumed  a constant  IJ. 
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Also,  boundary  conditions  for  transient  response  may  very 
well  change  with  time,  which  is  not  allowed  for  in  the  present 
formulation.  Such  problems  would  have  to  be  resolved  before  well 
founded  transient  results  could  be  calculated. 


CHAPTER  V 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Summary 

The  flow  in  a simplified  laminar  fluidic  amplifier  was  numer- 
ically and  experimentally  investigated.  The  study  divided  naturally 
into  three  major  parts. 

The  first  part  examined  the  numerical  simulation  of  a two 
dimensional  fluidic  amplifier  using  finite  difference  techniques  to 
predict  the  internal  flow  distributions  under  various  operating  condi- 
tions. Distributions  for  stream  function,  vorticity,  and  the  corre- 
sponding fluid  velocities  were  obtained  by  solution  of  finite 
difference  formulations  to  the  time  dependent  vorticity  transport 
equation  and  to  the  stream  function  equation. 

Appropriate  boundary  conditions  were  presented  and  examined 
for  a simplified  amplifier  geometry.  The  model  was  developed  such 
that  flows  with  different  Reynolds  numbers,  receiver  outflow  rates, 
and  deflections  of  the  inlet  flow  could  be  simulated.  A variety  of 
solutions  were  presented  as  contour  plots  of  the  stream  function  dis- 
tributions. 

The  predicted  flow  patterns,  and  their  changes  with  varying 
conditions,  are  reasonable.  The  numerical  routines  are  stable,  giv- 
ing convergence  for  all  cases  tried,  and  are  relatively  easy  to  use. 
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They  exhibit  no  dependence  on  the  initial  conditions  with  respect  to 
stability  but  do  have  time  step  limitations. 

In  the  geometry  chosen  the  flow  divides  into  a main  flow 
having  a varying  width  and  one  or  more  recirculating  vortices  which 
develop  on  each  side  of  the  main  flow.  For  the  lower  Reynolds  num- 
bers considered  the  main  flow  diffused  as  it  crossed  the  cavity  and 
relatively  weak  recirculation  vortices  were  formed.  At  the  higher 
Reynolds  numbers  the  main  flow  remained  in  a narrow  region  around 
the  centerline  of  the  amplifier.  In  each  side  of  the  amplifier 
cavity  this  narrow  central  flow  drove  one  main  reci rcul ation  vortex, 
which  in  turn  drove  one  or  more  smaller  vortices.  A characteristic 
of  the  flow  having  an  important  implication  on  the  pressure  solutions 
is  that  the  flow  becomes  more  complex  in  the  splitter  region,  i.e.  the 
gradients  of  the  flow  variables  become  larger,  as  the  deflection  of 
the  main  flow  increases  and/or  the  outflow  from  the  receivers  becomes 


more  unsymmetric. 

As  formulated  the  numerical  model  is  capable  of  predicting 
transient  flows.  However,  the  correct  interpretation  of  the  initial 
flow  field  as  an  initial  condition  for  a transient  flow  needs  further 
clarification.  Further,  the  numerical  simulation  does  not  predict  the 
"edgetone"  oscillations  which  were  experimentally  observed  at  higher 
Reynolds  numbers.  This  is  probably  due  to  inadequate  modeling  of 
the  boundary  conditions  and  the  use  of  a grid  size  too  coarse  to 
accurately  detect  small  movements  of  the  edge  of  the  jet. 
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« The  second  part  of  the  work  was  the  design,  construction,  and 

use  of  a large-scale  experimental  model  of  the  amplifier  considered 
in  the  numerical  simulation.  By  injecting  dye  into  water  moving 

B ‘ 

through  this  model  at  measured  flow  rates,  the  flow  patterns  were 

* made  visible  and  were  photographed.  This  was  done  for  flow  rates 

corresponding  to  Reynolds  numbers  of  20,  100,  and  290.  These  photo- 

I 1 

graphs  were  then  compared  with  the  appropriate  predicted  flow  fields. 

The  dye  patterns  illustrated  the  dividing  and  center  stream- 
lines  for  the  main  flow  and  the  flow  in  the  recirculation  regions. 

Dye  patterns  obtained  for  the  flow  having  a Reynolds  number  of  100 
were  particularly  clear  and  gave  good  agreement  when  compared  with 
the  predicted  flow  patterns. 

The  third  part  of  this  work  was  the  computation  of  pressure 
distributions  for  each  of  the  numerical  flow  field  solutions.  This  in- 
volved numerical  solution  of  a Poisson  type  partial  differential  equa- 
tion with  Neumann  boundary  conditions.  Problems  with  compatibility 
between  the  calculated  pressure  gradients  at  the  boundaries  and  the 
forcing  function  of  the  Poisson  type  equation  were  discovered.  This 
inconsistency  problem  was  examined  and  a method  was  presented  for  com- 
pensating for  the  inconsistency.  An  appropriate  convergence  criterion 
was  considered.  Converged  pressure  solutions  were  graphically  pre- 
sented as  plots  of  lines  of  constant  pressure.  Estimates  of  the 
accuracy  and  possible  shortcomings  of  the  predicted  pressures  were 
presented.  The  pressure  results  were  combined  with  the  flow  solu- 
tions to  generate  characteristic  design  curves. 
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B.  Conclusions 


Based  on  the  present  work  certain  conclusions  can  be  drawn 
with  regard  to  the  flow  field,  the  pressure  field,  and  their  simula- 
tion. These  are: 

1.  The  finite  difference  techniques  used  in  this  study  are 
a viable  means  for  modeling  the  steady  state  fluid  flow 
in  laminar  fluidic  devices  having  a relatively  simple 
geometric  design.  Therefore,  this  work  constitutes  a 
solid  basis  for  further  work  in  finite  difference  model- 
ing of  fluidic  amplifiers. 

2.  The  experimental  apparatus  built  in  this  study  allows 
visualization  of  the  flow  patterns  in  large  scale  ampli- 
fier models.  The  experimental  streamline  patterns 
qualitatively  confirm  the  numerical  streamline  predic- 
tions. 

3.  The  formulation  of  the  pressure  problem  as  a Poisson 
type  partial  differential  equation  with  Neumann  boundary 
conditions  is  insufficient  to  completely  specify  the 
problem  and  to  guarantee  adequate  numerical  modeling  of 
the  physical  process.  Also  needed,  and  often  neglected 
in  engineering  formulations  of  the  problem,  is  a state- 
ment of  the  consistency  requirement. 

4.  The  finite  differencing  necessary  to  obtain  the  boundary 
conditions  and  the  forcing  function  (s.,  of  Eq . (2.27)) 

' J 


cart  lead  quite  easily  to  violation  of  the  consistency 

requirement.  Although  the  differencing  is  done  correctly 

the  actual  magnitudes  produced  can  cause  the  violation 

due  to  truncation  error.  Further,  care  must  be  taken 

2 

with  the  differencing  formulas  selected  for  the  7 P 
terms  in  order  to  insure  consistency.  Consistency  of 
this  differencing  can  be  assured  by  the  methods  discussed 
in  Chapter  11,  Section  D.  These  methods  will  also  give 
the  actual  magnitude  of  any  inconsistency  due  to  trunca- 
tion error.  The  inconsistency  problem  can  be  allowed 
for  as  discussed  in  Chapter  IV,  Section  D.  The  resulting 
pressure  distributions  appear  acceptable  for  the  symmet- 
ric case  but  the  accuracy  seems  to  deteriorate  as  the 
flows  become  unsymmetric.  This  is  probably  due  to  resolu- 
tion problems  in  the  splitter  region  where  for  unsymmetric 
flows  the  gradients  become  steep  and  the  number  of  nodes 
is  small.  It  should  be  emphasized  that  use  of  the  alter- 
native line  integration  method  would  suffer  from  the  same 
resolution  problems  that  lead  to  inconsistency  in  the 
Poisson  pressure  formulation. 

5.  Problems  with  the  pressure  solution  produce  a restriction 
on  application,  at  this  point  in  time,  of  the  mathematical 
model  in  generating  amplifier  characteristic  curves. 
Therefore,  the  use  of  the  model  as  a design  tool  hinges 
on  future  refinement  of  the  pressure  solution. 


131 


6.  Assuming  the  pressure  problem  can  be  solved,  another 

consideration  in  use  of  the  model  is  the  cost  of  obtain- 
ing solutions.  The  computer  time  required  to  obtain  e 
set  of  solutions,  as  discussed  in  Chapter  IV,  Section  G, 
is  not  excessive,  but  it  is  not  trivial  either.  A 
characteristic  curve  was  generated  for  Reynolds  number 
of  100  in  6768  seconds  using  the  CDC  6600  system  at 
The  University  of  Texas  at  Austin. 


C.  Recommendations 


The  most  pressing  need  in  continuation  of  this  work  is  to 
solve  the  pressure  difficulties  encountered.  Based  on  the  work 
already  done,  it  appears  hat  insufficient  accuracy  in  the  differ- 
ence approximations  for  c regions  of  the  flow  field  is  the 

primary  cause.  Therefore,  the  following  recommendations  are  made. 

1.  The  next  step  in  this  work  should  be  an  effort  to  obtain 
higher  resolution  in  these  regions.  This  could  be  done 
by  either  the  use  of  a grid  with  variable  mesh  size  or 
by  use  of  some  type  of  grid  transformation.  Either 
approach  could  be  designed  to  give  a large  number  of 
nodes  in  the  regions  around  the  splitter  and  the  ampli- 
fier inlet  whore  gradients  are  large. 


2.  The  above  step  should  be  implemented  in  two  stages.  The 
first  should  be  examination  of  the  chosen  approach  on  a 
relatively  simple  flow  problem.  This  flow  should  not  be 
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that  of  a fluidic  il  ier  arid  should  have  a very  simple 
geometry.  Th<  piiv1'-  t.<  'o  would  be  to  reduce  the  varia- 
bles which  might  effect  the  pressure  and  flow  solutions. 
The  effects  of  the  grid  changes  on  the  pressure  solutions, 
particularly  ■ ■ jn  tude  of  the  inconsistency,  could 

then  be  more  easily  examined.  At  the  same  time  estimates 
could  be  obt  ined  f >i  what  effects  the  changes  would  pro- 
duce on  the  flow  ; . ■ d olution.  Naturally,  a flow  case 
for  which  a press  re  solution  is  available  for  comparison 
to  the  numerical  results  would  be  desirable.  The  second 
stage  would  be  impl  i he  refinement  on  the  full 
amplifier  model.  fully  the  programming  complexity 

would  be  reduced  by  the  experience  gained  in  the  first 
stage. 

3.  In  co i work  on  the  pressure  solution 

it  is  ecoi...  r : .0  i.,.c  pressure  values  be  obtained  experi- 

mentally foi  iiyi.tJ  the  numerical  predictions. 

4.  Once  the  prt  . *e  problem  s resolved  it  is  recommended 

that  the  po  apabilities  of  the  model  be  further 

explored.  The  abilities  to  vary  Reynolds  number,  outflow 
rates,  and  inlet  deflection  have  been  examined.  The  next 
step  would  he  to  examine  the  effects  of  varying  the 
geometry  of  the  amplifier.  This  miqht  include  varying 
the  receiver  widths,  the  receiver  angles  (presently  45°), 
or  the  ratio  ter  distance.  Finally,  the 
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geometry  ano  i .;.i . ; ,hould  be  added  to 

the  model. 

5.  It  is  recoi  he  model  to  predict 

transient  response  bo  t , lured.  In  particular,  interpre- 
tation or  in!  tic  , r.  ..  u oc  examined.  Also 

needed  is  the  de  elopmen  oi  boundary  conditions  which 
would  allow  edgetone  oscillations  to  exist  in  the  flow. 

6.  Finally,  it  is  suggested  that  boundary  conditions  and/or 

iterative  procedure  r.. eloped  which  would  allow  model- 

ing of  vented  a:,  lifur.-,.  As  discussed  in  Chapter  II 
such  flow  ant  flows  being  un- 

known in  either  dire  ' ;on  or  magnitude.  Unfortunately 

no  recommendat  ion  • ..  L made  as  i o who  c boundary  condi- 
tions might  be  used  to  a 4c  uately  model  the  vented 
ampl i fier. 
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APPENDIX 

Differencing  for  Boundary  Conditions  and  Forcing  Function  for  Pressure 
Solution 

As  discussed  in  Chapter  II,  in  the  calculations  for  the  pres- 

6P 

sure  distribution  values  must  be  obtained  for  — 

°y 


, SP 1 

and  -j — on 

Sx  i,j 


the  boundaries  and  for  s.  . at  all  nodes  in  the  region  of  interest. 

i ,J 

Boundary  gradients  are  obtained  through  appropriate  differencing 
of  Eq . (2.12).  Likewise,  s.  . values  are  obtained  by  differencing 

1 5 J 

of  the  relation  given  in  Eq . ( 2 . 29 ) . 

Selection  of  appropriate  difference  expressions  is  compli- 
cated by  the  fact  that  expressions  are  needed  for  wall  nodes,  where 
differencing  must  be  one  sided,  and  for  sharp  corners,  where  vortic- 
ity  values  are  discontinuous.  At  these  points  several  possible 
choices  for  difference  equations  were  available.  Various  combina- 
tions were  tried  in  an  effort  to  reduce  the  inconsistency  of  the 
pressure  formulation.  The  resulting  difference  equations  are  pre- 
sented in  Table  A. 1 . 
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Table  A.l.  Difference  Expressions  for  rP,  — , and  s.  . 

'■x  ny  1 , j 


Note:  Asterisk  denotes  multiplication 


2 < j < 40 


i = 11,  j = 41 


Right  boundary 
i = 31 , j = 1 


Location 

Difference  Equation 

Left  boundary 

i = 11,  j = 1 

5P 

s _L  * 

(~  '11,3  + 4 r’ll,2  " 3 ^llj) 

6x 

Re 

2 Ay 

SP 

■ =1* 

<‘r;12,l  ‘ ^10,1  ^ 

5y 

Re 

2Ax 

i = 11, 

6P 

: J_* 

(cn,j+i  " cn,.i-i) 

Sx 

Re 

2 Ay 

- u * (U12,J  " U1Q»J) 
11, j 


1 i+i  + i_i  " 2 'P 


- v 


ll.j 


n ,j+i  T "n  ,j-i  : ‘ "11, j 


Ay 


r 


SP 

s -1  * 

cll, 39  + 4 ?11,40  " 3 ^11. 41* 

Sx 

Re 

2Ay 

SP 

- ll  * 

(r12,41  " ’10,41 * 

<sy 

Re 

2Ax 

SP 

, _L* 

431  ,3  + 4 c31,2  " 3 ?31.1> 

Sx 

Re 

2 Ay 

SP 

=1  * 

(r32,l  * ^30,1 ^ 

5y 

Re 

2ax 

j 
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I 


L. 


Location 


Table  A. 2.  (continued) 

Difference  Equation 

■i  " Aoi  , _ i ) 


1 = 31, 

2 < j < 40 


<SP 


U 


i = 31 , j = 41 


6x 

Re 

2Ay 

' U31 

* (u32,j  ' 

U30,j) 

J 

2Ax 

* ^*31  ,j+l 

+ *31,3-1  ' 2 *11  .j* 

“ V31 

,j 

Ay2 

6P 

. -1  * 

(- 

C31 ,39  + 4 

^31 ,40  " 3 C31 ,41 ^ 

6x 

Re 

6P 

. -1  * 

(s 

32,41  ' ?30, 

41  ^ 

6y 

Re 

2Ax 

Lower  cavity  boundary 


12  < i < 18 

“ ' , j - 1 

24  < i < 30 

6P  _ 

-1  * 

^i+1,1  ' ci-l  ,1  ^ 

<5y 

Re 

2Ax 

i = 19,  j = 1 

6P 

-1  * 

C21  ,1  + 4 C20,l  ' 3 ^19,1 ^ 

6y 

Re 

2Ax 

i = 20,  j = 1 

6P 

-1  * 

^21,1  ‘ ?19,1* 

<5y 

Re 

2Ax 

+ u20 

* (*21 ,1  + *19,1  ‘ 2 *20,1 ^ 

+ w 

* ('  2ui9,l  " 3 u20,l  +6u21,1 

'20,1 


6Ax 
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138 
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Location 


Table  A. 2.  (continued) 

Difference  Equation 


6P 


1 * ^25 ,43  + 4 **25 ,42 


J "25,41 


5x 

Re 

2 Ay 

rB 

**25,41 

= - 

1 

2Ay2 

* 7 *25, 

41  + 8 *25 

1 

o 

•» 

*25 

6P 

i 

-1 

(- 

*“27 ,41  + 4 

^26 ,41  ' 

3 4 

,41 } 

<5y 

- 

Re 

2Ax 

6 P 

-1 

(- 

C25,39  + 4 

**25 ,40  * 

3 4 

,41 } 

6x 

+ 

Re 

2 Ay 

Left  receiver 

42  < j < 44 

w = i value  for 
wall  node 

Left  wall 


6P 

6x 


1 . (-S,.it2t4ew..Hl-3cW..1> 
Re  2Ay 


Right  wall 


6P 

<5y 


6P 

6x 


-1 

Re 


2Ax 


=1  * gw„j-2  + 4cw,j-l  ~ 3gw,j} 
Re  2 Ay 


6P 

<5y 


1 , (-V2..it4Vi..i-3^,j) 

Re  2 Ax 


45  < j < 56 
Left  wal  1 


6P  = _ <$P  = * (cw,j  ' VlJ+1^ 

6x  ’ Sy  Re  2Ax 


Right  wal 1 


j = 57,  i = 2 


j = 57, 

3 < i < 4 

j = 57,  i = 4 


Right  receiver 

42  < j < 44 

w = i value  for 
wall  node 

Left  wall 


Right  wall 
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L. 


Location 


Table  A. 2.  (continued) 

Difference  Equation 


fiP 

l 

* V2,j  +4  Vj  j ' 3 ’w. 

6y 

Re 

2Ax 

45  < j < 56 

Left  wall 

fiP 

fiP 

1 

* ^w,j  ' cw+l  ,j-l  ^ 

6x 

<5y 

Re 

2Ax 

Right  wall 

5P 

6P 

1 

* ^w-l,.i+l  " cw,j^ 

6x 

<$y 

Re 

2Ax 

j = 57,  i = 37 

6P 

fix 

37,57 

: §L 

fix 

36,56 

fiP 

fiP 

<5y 

37,57 

fiy 

36,56 

j = 57, 

fiPJ 

fiP 

38  < i < 39 

<5y 

i ,57 

fiy 

37,57 

j = 57,  i = 40 

fiP 

fix 

40,57 

fiP 

fix 

40,56 

fiPJ 

fiP 

5y 

40,57 

<Sy 

40,56 
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I 


Location 


For  all  wall  points 

20  < i < 22, 
j = 1 


11  < i < 31, 
2 < j < 40 


All  non-wall  points 
41  < j < 56 

Left  receiver 


Table  A. 2.  (continued) 


Difference  Equation 


■ 0 


S i » j 


(*i+i,i  * *i-i,i  ' 2 »i.i} 


2.0  * 

Ax" 


[“1+1.1  - ui-l  .l1 

2 

2Ax 

> 

Sf,j  * 2'°* 


(»1*1  ,.i  * Vl  ..1  - 2 *1  , j > 


AX 


Vu  ~ 


2Ax 


Si . j 


= 2.0 


^i-l,j-H  + *i+l,j-l  ' 2 ^i,j} 
2 


2Ax 


(<!>,• 


?1  J 


i-U+1  + ^i+l  ,j-l 


2Ax 


Table  A. 2.  (continued) 


Location 


Difference  Equation 


Right  receiver 


Left  receiver 
j = 57 

Right  receiver 
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